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LEE SIDE PLOW FOR SLENDER DELTA WINGS 
OP FINITE THICKNESS* 

Joachim SzodrUch 


Summary /1»* 

An experimental and theoretical Investigation was carried 
out to determine the lee side flow field over delta wings at 
supersonic speeds. The experiments were performed with models 
of the same slenderness s/1 =0.3 but different cross-section 
shape. The known types of flow, separated by the Stanbrook- 
Squire boundary into leading edge and shock-induced separation 
have been verified. However, further types of flow exist and a 
detailed survey of the boundaries is necessary. 

The influence of the cross-sectional shape on the lee side 
flow is discussed for the thick wing inside the region of shock- 
induced separation. In detail parameters like wedge angle at the 
center line, angle between upper and lower surface as well as the 
lower side shape are considered. 

A theoretical method to describe the flow field is lined out, 
where boundary conditions as a result of the experimental study 
are needed. The computed flow field with shock-induced separation 
is satisfactory. 


*ILR Report 23, Institute for Aviation and Spaceflight, Technical 
University, Berlin, 1977. ! 

##Numbers in the margin indicate pagination of original foreign 
text. 
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a angle of attack 

Y trace angle of Internal shock 

6* displacement boundary layer thickness 

e density ratio over shock 

\ wedge angle 

0 form angle (wedge angle for delta wing in the 

symmetry plane) 

k ratio of specific heats 
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t volume parameter 

$ inclination angle of shock 
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1. Introduction 

The aerodynamic design of hypersonic aircraft led to slender 
configurations with delta-shape plan form, strong sweepback of 
the leading edges and a substantial relative thickness, The aero- 
dynamic properties of such aircraft shapes are determined primarily 
in the case of supersonic flight by the flow conditions along the 
lower side. The lift is produced almost exclusively by the over- 
pressure along the underside. Suction forces along the upper side 
(lee side) have a very small influence. This makes it understand- 
able why only a small degree of attention was given to the flow 
field along the lee side of delta wings in the supersonic range. 
There were only a few surprising results in the determination of 
the aerodynamic heating on the top side of reentry bodies and 
space shuttle models which led to work on lee side flow. Measure- 
ments in the wind tunnel and in free flight tests [1-8] showed that 
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the temperature values on the upper side of a reentry body have 
a maximum value for zero Incidence angle (a « 0°) at the tip. 

Along the symmetry plane, the temperature decreases according 
to a hyperbolic law to the trailing edge. However, completely 
different conditions exist when the body has an Incidence angle. 

For example. If the Incidence angle Is a * 20°, there Is strong 
temperature drop Immediately behind the tip, which Is followed 
by a temperature Increase. The maximum value corresponding to 
a 8 0° is then exceeded. Within 20$ of the wing chord > tempera- 
ture differences of AT = 250° K have been measured. Depending 
on the Incident flow conditions, several temperature peaks could 
be determined. 

The results contrast to a certain extent with conventional 
ideas about the flow conditions along the lee side of delta- 
shape bodies. Up to the present It was always assumed, for the 
most part, that conical flow conditions prevail over delta wings, 
l.e., the state variables remain constant along rays through the 
tip of the wing. If one considers a continuous boundary layer 
development, one expects only slight and basically equivalent tem- 
perature variations. The results of the heating measurements show 
that the flow along the lee side of delta wings was important. 

Previously known investigations about the 1 ee side flow re- 
ferred primarily to thin delta wings . They showed that as the 
incident Mach number increases, there is a change In the type of 
flow — from a flow with a leading edge vortex to a flow with 
vortices caused by shock-induced boundary layer separations. This 
type of flow with shock-induced separations, which only occurs at 
the higher incident Mach numbers, is especially interesting because 
of the practical heating problems. 

The purpose of the present paper was to make a contribution to 
the problem of understanding the processes along the lee side of 
delta wings. The experiments attempt to make a close analysis about 
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the transition from the leading edge to the shock-induced separa- 
tion. By using thick delta wings with different cross-sections, 
we wish to investigate the Influence of different geometric 
parameters on the flow field. Using the experimental results, we 
developed a computation method with which one can describe the 
lee side flow for shock-induced separation. 


2. Present status of the results of lee side flow 


The flow field of a delta wing with an incidence angle in a 
supersonic flow can be broken down into two characteristic regions. 

The pressure side or underside faces the Incident flow and is 
dominated by the influence of the leading edge shock. The lee 
side or top side is in the wake of the model, and this is the 
region which will be described in the following according to our 
present state of knowledge. 

Figure 1 shows four flow models for the lee side of thin, 
slender delta wings with straight and sharp leading edges . The two 
figures on the left show the conditions for supersonic flow of the 
leading edges for both a small and a large angle of attack of the 
wing. The two sketches on the right show the corresponding condi- 
tions for supersonic incident flow of the leading edges. We can 
then give the following general description of the lee side flow 
according to Squire [9]. For Mach numbers normal to the leading 
edge which are substantially below the speed of sound and for small /9 
angles of attack, the flow is not separated along the entire top 
side (Fig. la). By increasing the angle of attack, leading edge 
flow separation occurs with conically, rolled-up vortices and 
secondary separations ('Fig. lb). 

Completely different conditions are found when the flow at 
the leading edge occurs at Mach numbers which are close to the 
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speed of sound or even above It. Figure lo shows that for small 
angles of attack there Is an expansion around the leading edge 
with nonseparated flow along the top side. A weak shock which 
runs through the tip of the delta wing Is responsible for the 
deflection of the flow. At higher angles of attack, the shock 
Is so strong that It Induces a separation of the boundary layer 
with conical vortices again (,Fig. Id). In both cases, the separa- 
tion in the central region of the wing shows a reattachment of the 
flow. 


Figure 2 shows the flow regions of thin delta wings as a func- 
Won of the angle of attach ana the Mach number M„, both normal 
to the leading edge. Four flow regions can be distinguished here, 
which were first discussed by Squire [10] in this diagram with the 
boundaries shown. Vortex formation at the leading edge and shock- 
induced separation with attached and separated leading edge shock 
are known phenomena. 

The Stanbrook-Squire region Indicates where the transition 
from the leading edge separation to shock-induced vortex formation 
occurs. This boundary was determined from experiments with models 
of different cross-section shapes and sweepback, but with a sharp 
or a round leading edge. For small angles of attack, the results 
agree with those of Squire [9] using the investigations of Lindsay 
and Landrum [11] for profiles with sharp or round leading edges 
and the Mach number range 0.6 < M^ < Q.8. The agreement is good. 

Considering the papers discussed above, we can see that the 
two flow types, leading edge separation and shock- induced separa- 
tion, are very important. Therefore, it seems appropriate to dis- 
cuss ’these in more detail in the following and to present existing 
theories. Furthermore, we would like to discuss a flow field which 
is very important for this work. At very high angles of attack, 
the vortex flow above the wing can explode and can produce large 
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area turbulent fields. This "vortex breakdown" is of Interest 
especially because of certain experimental results In this paper. /10 
In the last section we will briefly describe aerodynamic aspects 
In the heating region along the lee side, because Important In- 
formation about the top side flow will be given from It. 

2.1. Leading edge separation 

The leading edge separation of delta wings with a sharp lead- 
ing edge can be explained In the following physical terms . There 
Is a flow around the leading edge from the pressure side to the 
leeward side. In this way, a high negative pressure forms in the 
region directly adjacent to the leading edge. The large pressure 
gradient which Is produced in this way between the top side and 
the bottom side produces a flow separation. It starts from the 
leading edge and develops Into a spiral conical vortex above the 
lee side. A flow component Is produced along the bottom side of 
the primary vortex, which Is aligned with the leading edge. How- 
ever, It is not large enough to overcome the pressure increase be- 
tween the underpressure peak Induced by the vortex and the leading 
edge. There Is a secondary separation which occurs, and again it 
rolls up Into a conical vortex. Figure 3 shows the entire flow 
field together with the characteristic stream lines and the pressure 
distribution. 

First, we will discuss the vortex as the main element of the 
lee side flow of leading edge separation. According to measurements 
of Earn shaw [12], the primary vortex can be primarily separated into 
three parts. First of all, there Is a vortex surface which emerges 
from the leading edge and which "supplies" the second part, the vor- 
tex Itself. This can be assumed to be frictionless and conical. 

The third part Is the vortex chord, a small region in the center of 
the vortex In which friction influences are not negligible. In this 
region, large gradients of total pressure, static pressure and 
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velocities have been measured. Experiments show, In comparison 
with the theoretical analysis of the vortices by Hall [13] » that 
there Is good agreement with respect to the vortex shape and 
velocity distribution. 

In general, this vortex structure was confirmed In most ex- /II 
perlments, for example, the flow field Investigations of Drougge 
and Larson [It], as well as Thomann [15]. One detail effect we 
wish to Investigate Is how the vortex position and pressure dis- 
tribution or life coefficient change under the Influence of vari- 
ous parameters. 

According to Okerbloom and Sarantsev [16] , an Increase in 
the Mach number for fixed angle of attack leads to a displacement 
of the vortex centers from the leading edge towards the central 
plane. On the other hand, the height of the vortex core above 
the model surface remains about the same. Prom this behavior, 
we can see that the underpressure peak also decreases with In- 
creasing Mach number and angle of attack. Since the size and 
position of the underpressure maximum varies greatly In the vari- 
ous experiments and theories. Hurley [17] attempted to find a 
correlation of the various data, using the angle of attack, sweep 
angle and Mach number. 

Influences of geometry of delta wings on the flow conditions 
have been reported by Fellows and Carter [18] for wings and wing- 
body combinations. Greenwood [19] showed that with differing 
sweep angle, large negative pressure coefficients are Induced near 
the leading edge for all Mach numbers, and there is only a slight 
change In the central region of the delta wing. 

According to Fig. 1 and oil film image photographs, there 
are wall stream lines parallel to the incident flow in the region 
of the attach flow near the central plane. In contrast to this, 

Morris and Couch [20 ] found In wind tunnel tests an additional 
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component In the direction of the symmetry line in the central 
region after reattachment of the flow with very thick reentry 
bodies. The oil lines then run tangentially to the central 
plane downstream. 

4 

In the case of asymmetric vortex separations from the top 
side, one finds a further deviation from the normal behavior of 
the flow, which can occur primarily for wings in a symmetric. 

Incident and slipping flow. Gapcynskl [21] investigated this 
with elliptical cylinders, for which the asymmetric states were 
determined when a certain angle of attack was exceeded. /12 

In general, most results show a high similarity of the flow 
conditions in the subsonic and supersonic range when there is 
leading edge separation, as experiments of Lee [22] have con- 
firmed. 

Most of the calculation methods for the slender delta wing 
with subsonic leading edg are based on the theory of Jones [23]. 

The flat delta wing is treated in planes perpendicular to the 
incident flow, just like a plate with vertical Incident flow, 
using potential theory. In this first paper, the leading edge vor- 
tex was not included; therefore, the following theoretical models 
of Legendre [2*1], Adams [25 J and Edwards [26] were concerned with 
the extension of the problem to two concentric vortex cores above 
the lee side in a potential flow. Brown and Michael [273 intro- 
duced a further substantial improvement of the vortex model, by 
connecting the concentric vortex cores with the leading edge using 
a vortex surface, and the vortex cores are represented as vortex 
lines starting from the tip of the model. Mangier and Smith [28] 
extended this model further. In this case a spiral-shaped vortex 
is divided into an outer part which emanates from the leading edge 
and an inner part which forms the foundation of the calculations. 
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The two last models, however, did not bring about any substan- 
tial Improvement In spite of the fact that different models were 
used, so that both were used as the foundations for further cal- 
culations of Smith [293 and Nennl and Chee Tung [30]. The formed 
and thin delta wings were calculated by Carafoll and Stalcu 1 3ll , 
by superimposing known flow fields and flow fields calculated 
with the method shown above. A refinement of the vortex model Is 
achieved by simulation of the secondary vortex, which Is re- 
placed by a stagnation point flow by Pershing [ 32 ] . 

The above models for calculating the leading edge separation 
are based on the theory of slender bodies and Mach number influ- 
ences are not considered. In experiments it was shown that as the 
Mach number increases, the vortex becomes flatter and Its height 
over the leeward side decreases. These results were used by 
Ktiehemann [ 331 and Squire [3*U and they replaced the real flow /13 

by vortices which are directly above the top side of the wing In 
their theoretical calculations. 

Polhamus [35] gives a completely different solution method. 

The method is based on the "intuitive" assumption that the normal 
forces in the case of leading edge separation with flow reattach- 
ment in the central range are equal to the suction forces when 
there is flow around the leading edge with reattachment on the 
top side. This trial solution also considers Mach number Influ- 
ences. Very similar equations are obtained for the lift coeffi- 
cient from a mathematical model of Coe [36], even though here the 
model of concentrated vortex cores is used again. Along the axis 
of these vortex lines, the effects of mass supply are introduced 
in the calculation using additional sinks. Even though the results 
of Polhamus agree much better with experiments compared with those 
calculated by Coe, the latter method seems to justify the "intuitive" 
procedure of Polhamus. An additional method of calculating conical 
flows is the method given by Hummel [373 , which determines the flow 
variables according to the theory of slender bodies in cross-sec- 
tional planes with source distributions and vortex distributions 
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along the contour. In this method, arbitrary body shapes can be 
used, but the leading edge vortices were not considered in the 
flow model of this theory. 

Table la Is a summary of the various theories for leading 
edge separation which deals with new flow models. 

Summarizing, we can say that there Is a vast amount of 
knowledge about lee side flow in the case of leading edge separa- 
tion. In experiments, it was possible to clearly Identify the 
dominating vortex systems. However, simplified flow models were 
introduced in the theory. The agreement between the theory and 
the experiments regarding the pressure distribution and position 
of the vortices is partly unsatisfactory. 


2 . 2 . Shock-Induced separation 

Shock-induced separation along the top side of delta wings 
occurs in most cases of a supersonic leading edge. If the lead- 
ing edge shock Is attached, then it is possible to consider the 
top side and the bottom side separately. There is no flow around 
the leading edge, but Instead there is an expansion which can be 
well-described using the Prandtl-Meyer flow. In this region the 
flow is attached on the top side. The expansion then causes a 
deflection in the direction of the central line. Because of 
symmetry, the flow must be parallel to it in the central plane, 
l.e., the transverse speed must be zero there. Since this back- 
ward deflection is to occur in a supersonic flow, it Is created 
by a compression shock which runs through the tip of the delta 
wing. If the intensity of the shock is large enough, then there 
is a separation of the boundary layer downstream and a vortex 
forms which has a sti’ucture similar to the case of leading edge 
separation (Fig. 4). This model is supported by the investiga- 
tions of the flow field with special regard for the position and 
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appearance of the Internal shock, which was carried out by Bannlnk 
and Nebbeling [38-39]. The pressure measurements of Larcombe [40] 
clearly show the existence of the expansion with the subsequent In- 
ternal shock In the vicinity of the symmetry plane. Also, the 
pressure distributions of Squire [4i~k2] on delta wings with 
different cross-section shape and sweepback show this. Here, it 
was possible to show that changes in the cross-section geometry 
of the delta wings does not bring about any substantial changes In 
the overall flow field, which has also been confirmed by Investiga- 
tions on the space shuttle 18, 43-44]. 

Pike [453, in experiments with round and sharp leading edges, 
showed that blunt edges change the lift only slightly, but the 
high pressure in these regions does contribute to the drag. The 
lowering of the peak of the delta wing [10-12] to reduce the heat- 
ing also influences the flow in the symmetry plane. 

Rao and Whitehead [5] established a flow model for the better 
understanding of the processes at high Mach numbers. It is shown 
in Pig. 5 and shows that the flow Is not conical over a large part 
of the delta wing. The basis of this Is a hypersonic "two-layer" /15 

boundary layer with various momenta. The pressure increase through 
the internal shock is different for the two layers and in this way 
a vortex is created inside the boundary layer. The vortex again 
in turn brings about a rarefaction of the boundary layer in the 
central region, which is related to a reattachment of the flow in 
the inner part of the wing. The flow model of Rein [46] , shown in 
Pig. 6, is also not conical. It is the result of oil film photo- 
graphs. If the leading edge shock is separated, one can observe 
leading edge vortices in the tip region, which, however, have to 
be explained here as the result of the mutual influencing of the 
shock and the separation bubble. The separation bubbles then form 
downstream at the leading edge, and again internal shocks can occur 
on them. Since the vortices can no longer be "supplied" by the 
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leading edge, they appear as free vortices in the central region 
of the wing. Experiments with two-dimensional wedges have shown 
that at least qualitatively there is a similar behavior between 
the flows perpendicular to the leading edge of a delta wing and 
the flow around a wedge with the same aperture angle, at least 
in the region of the leading edge. 

Large wedge angles between the top side and the bottom side 
perpendicular to the leading edge or high angles of attack lead 
to a separated leading edge shock. Because of the subsonic region 
which is downstream from the shock, leading edge separation can 
again occur, according to Ghorai C 473 . Since the subsonic region 
is limited to a relatively small region in the vicinity of the 
leading edge, a shock-induced separation occurs, as was confirmed 
by investigations of Collis (48) with six delta wings. Figure 7 
gives a comparison of the flow models and shows the cross-section 
flow over a thick (Collis) and a thin (.Cross C 493 ) delta wing in 
the hypersonic range. The boundary layer shape is different and 
this is related to the number of internal shocks, which can be one 
or two along each half of the wing. In both models, the shape of 
the boundary layer in the region of the symmetry plane seems un- 
realistic, because here one can assume a continuous transition from 
one half of a wing to the other. The most important data of the 
experiments described above are collected in Table 2. 

In the following we will now present the most important 
theories for describing the lee side flow for supersonic leading 
edges- and will compare them with experiments. Maslen [50] dis- 
cusses the top side flow only as a special case of the general flow 
around a delta wing with the solution of the nonlinear equations 
for conical flow. The calculation has no shocks and is divided in- 
to Prandtl-Meyer expansion and hyperbolic or elliptic regions, 
within which it is possible to give a solution of the differential 
equation using characteristics or relaxation. The conical flow 
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over a plane delta wing with supersonic leading edge was solved 
by Clarke and Wallace [51] using an Integral method In second 
order terms with respect to the angle <>f attack. Internal shocks 
cannot occur here. On the other hand, Fowell 1 52] calculated both 
the continuous case without a shock as well as the discontinuous 
case which, under some conditions, is the only possible solution. 
Babayev [533 has proven that continuous flow cannot occur. This 
solution consists of Prandtl-Meyer expansion, determination of 
the shape of the internal shock and the shape of the adjacent 
central region using iteration methods. Pike I 543 calculated the 
pressure distribution along the top side of a surfboard using the 
linearized theory. Prom parameter studies, he concludes that with 
increasing sweepback of the leading edge, the internal shocks in- 
crease in intensity. In order to give a correct interpretation of 
this result, we have to mention again that the internal shock de- 
pends on the additionally induced component during expansion around 
the leading edge. The Prandtl-Meyer expansion, in turn, depends 
on the angle of attack and Mach number perpendicular to the lead- 
ing edge. For a large sweepback, however, the M N is reduced, but 
the angle of attack becomes larger, and in this way the estima- 
tion of the intensity of the internal shock is made more difficult . 

A three-dimensional characteristic method in which the internal 
shock is represented by isentroplc compression was established by 
Beeman and Powers [55] . The special difference method of Kutler 
and Lomax 156] is used in order to solve the three-dimensional, 
frictionless and nonlinear basic equations for the top side as well. 
An extension by Walkden, et al. 1573, using boundary values from 
the characteristic method, did result in improvements to the calcul- 
ation program, but the internal shocks with a real intensity are 
not calculated. Miyazawa 158], using comparisons between experi- 
ments and modified theory of Kutler and Lomax, showed that on the 
lee side of cones, no good agreement could be expected. In par- 
ticular, for the hyperbolic region of the cross-section plane. 
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Christophel [ 59] determines the flow field using the characteristic 
method. This seems to be limited to the localization of the In- 
ternal shock. To a limited extent, the "method of lines" of South 
and Klunker T60] can also be used for the lee side, but poor com- 
parison possibilities result for high angles of attack, because In 
the real flow friction effects prevail. Extensive flow field In- 
vestigations of Gross [493 were the foundation of a calculation 
method which, however, Is based In part on empirical data. A model 
Is used In which the flow expands at the boundary layer and it Is 
deflected by a shock. The shock which Is attached to the boundary 
layer produces a separation directly on the surface, which again 
results In a substantial thickening of the boundary layer behind 
the shock. The calculation method, however, is restricted only to 
expansion over the boundary layer, the internal shock and the ad- 
jacent central region. The vortex formation Inside the boundary 
layer is Ignored. The shape of the boundary layer and a parameter 
of the area equation of the Internal shock must be Introduced Into 
the calculation as empirical data. 

Table lb shows the various theories of shock- Induced separa- 
tion as well as the different flow models. 

Numerous experiments in the area of shock-induced separation 
have been carried out, but in only a few papers was it possible to 
detect the internal shocks. This already indicates that the knowl- 
edge about the processes in this type of flow is still deficient 
and is usually limited to high Mach numbers and thin and flat delta 
wings. The theoretical models which attempt to calculate the flow 
field with closed solutions have no shocks. Also, the few papers 
which include the internal shocks are not satisfactory compared 
with experiments. As far as is known, the model of Cross is the 
closest to the real flow, which was selected in the present paper 
as the basis of the theoretical analysis. 
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2.3. Transition from leading edge separation to shock-induced 
separation 
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Up to the present there Is no systematic Investigation about 
the transition from vortex formation at the leading edge to the 
shock-induced separation. This change In the flow type Is very 
likely continuous and as a first approximation can be assumed to 
occur within the Stanbrook-Squire region (Pig. 2). This region 
was determined by considering the different vortex formations from 
all of the available experiments, both directly at the leading edge 
as well as inside through an Internal shock. However, we should 
point out that It is difficult to carry out an exact analysis of 
the flow type because of ambiguities in the experimental Investiga- 
tion methods. For example, oil film photographs have very similar 
characteristics for leading edge separation and shock-induced 
separation. This could be the reason why a few experiments, for 
example, those of Whitehead and Keyes [1] , cannot be classified 
In the boundary region. 

Figure 8 gives a possible model on how the transition of the 
flow fields could occur. It was given by Squire [6l], supported 
by "vapor screen" photographs. Starting with shock-induced vortex 
formation, the flow perpendicular to the leading edge at small Mach 
numbers- is no longer capable of completely expanding. A separation 
at the leading edge occurs and, in addition, there is an expansion 
with subsequent internal shock. If the Mach number drops again, 
then the leading edge vortex is enlarged, whereas the expansion 
and the shock decrease in intensity, until finally only leading 
edge separation occurs with the secondary vortex. 

An additional representation of the processes at the bottom 
side which can be related to the change in the flow type along the 
top side was also derived by Squire [ 6l] from the "thin shock 
layer" theory. For a separated leading edge shock, the theoretical 
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analysis results In a discontinuous transition from a stagnation 
line In the symmetry plane to an additional one In the vicinity 
of the leading edge. The number of stagnation lines depends on 
the flow region, so that according to the leading edge (1 stagna- 
tion line) or the shock- Induced separation (2 stagnation lines), 
we can Identify the flow range. The projected stream lines of 
the bottom side are also shown In Pig. 8, and a continuous transi- 
tion was assumed there. 

The change from leading edge-induced separation to shock- 
induced separation and the related flow processes at the bottom 
side are only to be Interpreted as a working hypothesis, because 
no experimental investigations have been concerned specifically 
with this problem. 


2,4. Vortex breakdown 

The flow field on the lee side of a delta wing with leading 
edge separation can collapse for certain changes in the parameters, 
for example, if one Increases the angle of attack. This phenomenon, 
also called vortex explosion, was explained by Werle [62] as an 
expansion of the free, spiral-shaped vortex because of transition 
from laminar to turbulent flow. Brooke -Benjamin [63] also does 
not Interpret the collapse as an instability, as was stated In 
other investigations, but as a second dynamic state of the flow. 

The "vortex breakdown" usually starts downstream from the trailing 
edge and wanders upstream with increasing angle of attack. If a 
position Is then reached in the vicinity of the trailing edge, 
then even small angular changes are sufficient to influence the 
position of the bursting point. If the trailing edge Is exceeded 
and if the angle of attack is increased further, then according to 
Lawson [64], a stable position in the vicinity of the tip of the 
delta wing occurs and for sweep angles between A -70° and 85°, 
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they are about at 37 % of wing length. Peckham [65] determined 
that the bursting position depends on a combination of angle of 
attack and sweep angle. Since Inverse proportionality exists be- 
tween the two. In practice this means that "vortex breakdown" Is 
an upper limit with respect to angle of attack, and a lower limit 
with respect to sweep angle for the flow with leading edge separa- 
tion. 


The results of Wentz [66] , Sarpkaya [671 are of interest for 
the present experiments. In the critical range, the bursting of 
the vortices apparently Is very sensitive to disturbances in symme- 
try, whereas an Increase In the pressure in the direction of the 
vortex axis has a destabilizing effect. Even though relatively 
little is known about the process of bursting and only a few 
parameters can be recognized with certainty, Brooke-Benj amin [68] 
and Jones [ 69 ] performed a theoretical analysis of this problem. /20 


These phenomena are observed in all subsonic incident flows 
and only Elle [70] and Lambourne and Bryer [71] carried out ex- 
periments near the speed of sound. For supersonic Mach numbers, 
in [72-73] similar flow conditions were found in experiments. 

Here again the "vortex breakdown" was noticed because of pressure 
fluctuations and a decrease in the underpressure peak downstream 
of the bursting point. Oil film photographs of burst vortices in 
the subsonic and supersonic range have the same characteristics. 
From this we can conclude that supersonic flow shows essentially 
the same parameters and flow processes when there is collapse of 
the vortex. 

Here we have attempted to discuss briefly the special state 
of the flow over the delta wing, which can only occur for certain 
incident flow conditions. Even in the framework of the present 
investigation, the occurrence of the "vortex breakdown" is think- 
able, so that experimental results also have to be evaluated from 
this point of view. 
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2.5. Heating 


Already In the Introduction we mentioned that many problems 
of lee side flow are the result of heating Investigations. The 
high local temperatures which occur In the central region of the 
wing are of special Interest. Experimental Investigations of 
Whitehead and Keyes [1] led to a flow model for explaining these 
heating peaks. On the lee side, spiral-shaped vortices are 
created, and their circulating motions Induce components directed 
downwards In the vicinity of the symmetry plane. On the other 
hand, near the body surface, a tangential velocity component Is 
produced by the vortices which removes layers with low energy 
from the central region In the direction of the leading edge. 

The temperature peaks then result from the type of stagnation 
point flow In the symmetry plane and the boundary layer thickness 
which Is reduced there. Whitehead 12] was essentially able to 
confirm these tests and even found two peaks In the temperature 
distribution along the central plane. One maximum was In the 
vicinity of the vortex origin and the other was produced by transi- 
tion of the boundary layer. /21 

The complete suppression, or at least reduction, of the high 
local temperature peaks by forming of the lee side or the leading 
edge or by dropping the model tip was not completely possible 

13-5 3. 

The essential parameters which influence the temperature peaks 
were investigated by Whitehead, Hefner and Rao [6] and were summar- 
ized by Hefner T8] . In the lower hypersonic range (M < 6), the 
temperature maxima depend greatly on the Reynolds number, but Its 
local position can only be Influenced by the angle of attack. In 
addition, there is a limiting Reynolds number below which the heat- 
ing peak is reduced greatly and abruptly with decreasing Re number. 
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All of the previous experiments have shown that the greatest 
temperature Increases occur along the top side In the low hyper- 
sonic range. At higher Mach numbers, the reduction In the bound- 
ary layer thickness discussed above In the central region does not 
occur, and Instead there Is a thickening. In this way, the flow 
conditions change so that the temperature peaks are greatly sup- 
pressed. 

This summary does not assume to be complete, instead we only 
wanted to indicate the aerodynamic aspects In the area of heating 
problems . 


3- Experimental facilities 

The analysis of the flow field on the lee side Is essentially 
based on experiments in the wind tunnel using models with various 
design parameters . In the following we will present the wind 
tunnels, models and various test methods. 


3.1. Wind tunnels 


Three different wind tunnels were used for the tests discussed 
In this report. The first experiments were performed In the super- 
sonic wind tunnel of the Cambridge University Engineering Department 
(CUED). The tunnel is a blow-down type with a test section cross- 
section of 0.12 x 0.18 m and was operated at a Mach number M^ = 3-5 
and a Reynolds number Re * 5.0 x 10^ 1/m. 

Additional Investigations were performed In the supersonic wind 
tunnel of the Royal Aircraft Establishment (RAE) in Bedford with a 
0.9 x 1.2 m test section and for Mach numbers of = 3.5 - 4.5 


an <3 Reynolds numbers of Re * 2.2 x 10 ^ - 2.6 x 10*^ 1/m. 

The high speed wind tunnel of the Institute for Aerodynamics 
and Spaceflight (ILR) with a test section of 0.15 x 0.15 m was 
used in the speed range of M^ * 2.0 to 4.0 and for Reynolds num- 
bers around Re * 1,0 x 10*^ 1/m to test new and geometrically 
similar models. 

Since we expected large Interferences between the quasi- 
frictionless flow regions and the regions with friction. It Is 
appropriate to compare the order of magnitude of the Reynolds 
numbers for free flight and In wind tunnels. Figure 9 shows 
the Reynolds numbers of the wind tunnel models using the present 
report as well as those of SST aircraft with design Mach numbers 
between M^ ■ 2.2 and 3.5 as well as for two reentry bodies. The 
differences which occur are considered In the discussion of the 
experimental results. 


3.2. Models 


All of the models investigated are delta wings with straight 
and sharp leading edges and conical surfaces. Figure 10 shows 
the different cross-section shapes of the Models I. These bodies 
were measured In the Cambridge and RAE wln<2 tunnels and have a 
slenderness ratio of s/1 = 0.31. Therefore, five different top 
side shapes are available: two delta shapes and two conical 

shapes, as well as a flat leeward side. The wind tunnel Models II 
in Fig. 11 have a slenderness ratio of s/1 = 0.30 and were tested 
in the ILR wind tunnel. Two cross-section shapes are the same as 
for Model I, and in addition several experiments were carried 
out using the Nonweller wave rider. 

The slightly modified conical model shown In Fig. 12 has a 
shape which is probably close to a practical shape and it also has 
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a conical surface. 


The pressure traps In the span direction were located at a /23 
relative wing chord of x/1 ■ 0.55, U.7Q and 0.80. Further de- 
tails are shown In the figures. 

In all cases of this Investigation, the angle of attack Is 
measured with inspect to the ridge line or the highest line of 
the top side. 


3.3. Model suspension 

The type of model suspension In the wind tunnel can create 
a disturbance to the flow around the delta wing under some condi- 
tions. In order to estimate possible influences, we show the sus- 
pensions In the three wind tunnels In Fig. 13 • In the CUED tunnel, 
the wedge-shaped sting only penetrates Into the wake on the pres- 
sure side of the model. In the RAE wind tunnel, there Is a rela- 
tively thick suspension wedge, both on the pressure side and on 
the lee side. In the ILR tunnel, the relatively thin sting Is 
located about three model lengths x K /l = 3.0 downstream of the 
trailing edge. 


3-^. Test methods and errors 

In addition to static pressure measurements In the longi- 
tudinal and span directions on the model surface, Pitot pressure 
measurements were carried out at various heights above the body 
using a probe. It is difficult to estimate the influence on the 
flow field and the resulting measurement errors If the probe enters 
the region of the boundary layer or into a region of local subsonic 
flow. If one uses Pitot rakes, as in the present case, then we 
also have to consider the Interference between the individual 
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measurement tubes. Quincey and Calllnon [741, at a Mach number of 
M^ » 1.6, carried out experiments with rakes of different geome- 
tries. According to this, the Intermediate space should be at 
least three times the external diameter of the tubes in order to 
avoid Influences. The present experiments were performed with a 
rake where the spacing between the individual Pitot tubes was four 
times the external diameter, and so that there should be no inter- 
ferences even for local Mach numbers smaller than M e * 1 . 6 . 

In addition, the error which is the result of oblique incident /24 
flow to the Pitot tubes has to be estimated, because the probe had 
the same angle at all points of the flow field with respect to the 
incident flow. Figure 14 shows the calibration of a Pitot tube for 
two Mach numbers for angles of attack of up to « * 23 °. The Pitot 
pressure measured at different angles of attack deviates by a maxi- 
mum of 3$ compared with the value for a ■ 0. For small angles of 
attack of up to about a = 4 °, no Pitot pressure losses could be de- 
tected. 


In addition, errors occur in the measurements because of the 
nonuniform flow in the wind tunnel and also because of the measure- 
ment apparatus, such as the manometer, pressure transducer, scanner 
and pressure lines. Also, errors caused by asymmetry in manufactur 
ing the model or nonexact alignment in the wind tunnel can lead to 
different results on the two wing halves. 

An error estimation for the pressure coefficient in the CUED 
tunnel resulted in Ac = -J-.002 and even smaller values for the RAE 

p - 

wind tunnel. 

At low pressures, the measurement is much more problematical, 
but in the atmospheric ILR tunnel for static pressure of the in- 
cident flow of p = 34 torr, we only detected a maximum error of 

UU 

Ac p =+ ,0015 . If we consider all of the error sources, it seems 

that a deviation of a maximum of Ac = + .0025 seems realistic in 

P ~ 
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all tunnels . 


Investigation methods for the visualization of the flow in- 
cluded oil film photographs, shadowgraph photographs and Schlleren 
photographs . In the Cambridge wind tunnel we took oil film photo- 
graphs after each test run, but before we observed to make sure 
that no lai*ge changes had occurred b 6 cause the tunnel was turned 
off. In the ILR wind tunnel, all of the oil film Images were 
taken during the test even though here there Is the disadvantage 
of a reduced contrast and sharpness. Instead, In this method one 
is able to better recognise the dead-water regions in which the oil 
collects to form drops In some areas. 

Even though we cannot make any definite statements from the 
analysis of the oil film photographs, we can still make a compari- 
son of several characteristics. If we consider the vortex forma- 
tion at the leading edge with the secondary vortices and attached 
flow in the central range, then we give the definitions from an 
oil film photograph of Pig. 15. The positions shown have been 
examined using other investigation methods. When specifying the 
vortex Intensity, we assumed that the Increase in the vortex in- 
tensity was related to a large radial speed, which again makes 
the wall stream lines run at a larger angle <j!. 


4 . Discussion of experimental results 

Static and Pitot pressure measurements, as well as flow 
visualization methods, were used as the most important test 
methods in the experiments. The results will be used In the 
following to analyze the flow regions for the delta wings under 
discussion and to develop models of the processes on the lee side. 
We assume the already described two flow models for the leading 
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edge separation and the shock-induced separation. 


When analyzing the test results, it has been found to be 
appropriate to use the components of the Incident Mach number 
^ao and ^e angle of attack a in the plane perpendicular to the 
leading edge of the wings Instead of the values themselves. 
These components are found from 
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The angle 0 Is In the symmetry plane between the plane of 
the leading edges and the ridge line, or the highest line on the 
top side through the tip of the model. Therefore, for the flat /2 6 

lee side we have 0=0° and the equations given reduce to the 
equations for the flat delta wing. The equations are given In 
Fig. 1 6. The upper part of the figure shows the normal angle of 
attack for different top sides and the lower diagram shows the 
normal components of Mach number for the flat top sides as a func- 
tion of angle of attack. Figure 17 shows the regions for the 
angles ot^ and the Mach number M N which were selected In the ex- 
periments. The tests with delta wing for a flat top side coincide 
directly with the Stanbroolc-Squlre region; therefore, they can be 
subjected to a detailed analysis. For all other models, the 
and M n values are to the right of the Stanbroolc-Squlre region. The 
experiments with these models attempt to determine the influence 
of the cross-section area for shock-induced separation. Before we 
can discuss these problems, however, we must clarify whether conical 
flow can be assumed under the prevailing conditions. 
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4.1. Conical and non-conlcal flows 


Prom the present state of knowledge of the lee side flow, we 
can see that conical and also non-conlcal flow fields can be dis- 
tinguished. For most Investigations with conically shaped sur- 
faces, conical flow can be assumed as a good approximation. How- 
ever, Cross [49], Squire [4l] , based on the models of Whitehead, 
et al. [63 and Rein [463, assumed non-conlcal flow over most of 
the wing. In the following we will examine to what extent these 
statements apply for the delta wings and the Incident flow conditions 
Investigated here. ■ 

Figure 18 shows the result of static pressure measurements 
for six different wing chords on the delta-shaped lee side of 
Model I. The results are the isobars for two angles of attack. 

These figures will be associated with the corresponding oil film 

photographs and can be considered representative for all of the 

other models. For both angles of attack, there is approximately 

conical flow, but we have to point out that this assumption is 

only valid to a limited extent. Deviations from conical flow 

were found for very small angles of attack as well as for very 

large angles of attack, especially in the tip region or in the 

region of the trailing edge. Since the models were relatively /27 

small, it is not possible to exactly examine the tip region. 

Therefore, the investigation was concentrated on the central 
region, in which to a good degree of approximation there is conical 
flow for most of the incident flow conditions. It is assumed that 
in spite of this restriction, the results are representative for 
the entire flow field. 

For wind tunnel tests with fixed model suspensions, one has 
to especially consider the influence of disturbances in the wake, 
which, just like the tail flow, can induce a deviation from conical 
conditions in the region of the trailing edge. This was established 
in tests and is discussed in a special chapter. 
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4 .2 . Results on the transition from leading edge-induced 
separation to shock-induced separation 

The processes Involved in the transition from the leading 
edge-induced to shock- induced separation are, for the most part, 
unknown and will be discussed In the following. Based on an 
evaluation of experimental results. Squire [10] gave a transi- 
tion region called the Stanbrook-Squire region. Figure 17 shows 
that the models used here with the flat top side cover the 
Stanbrook-Squire region well for Mach numbers between M * 2.5 
and 3.5. The flat lee side was selected so that the normal Mach 
number is clearly in the supersonic region and also because there 
is sufficient comparison material from other investigations. 

First of all, we analyzed the flat top side of the delta model 
and its flow field on both sides of the Stanbrook-Squire region, 
in order to then describe the change in the flow types with in- 
creasing incident Mach number for a representative angle of attack 
of the wing. 

4.2.1. Separation with leading edge vortex 

Leading edge vortices, according to Fig. 17, can be expected 
at a Mach number of M QO = 2.5 and angles of attack between ar.= 5° 
and 15°. Very large angles of attack were first not observed be- 
cause the Stanbrook-Squire limit up to the present was only limited 
to values smaller than a = 50°. /28 

According to known experiments , the following model prevails 
for the lee side flow for leading edge separation. Directly at 
the leading edge, a primary vortex separates, which again induces 
a secondary vortex along the outer wing region. Near the symmetry 
plane, there is a similar flow which reaches this region from 
above the vortex. Because of these conditions, a constant pressure 
is established in the central region, whereas the primary and 
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secondary vortices produce underpressure peaks. This model was 
confirmed by our own experiments, which Is shown In the following 
for a delta wing with a flat lee side at M = 2.5 and a = 10°. 

CO 

Figure 19 shows the different available Information for this 
incident flow case. A direct comparison allows one to compare 
the statements many times. Attach flow In the central region was 
confirmed by oil film tests, Pitot tests and static pressure tests. 

The position of the primary vortex Is found at about y/s =0.55 
and z/s =0.20 and can be localized just like the secondary vortex 
at about y/s =0.78 and z/s = 0.08, as well as from oil film photo- 
graphs and from the different Schlieren photographs. Based on the 
low intensity of the secondary vortex and the closeness to vacuum 
conditions, one does not observe any second underpressure peak in 
the pressure distribution. Therefore, we can establish that the 
well-known ideas about leading edge separation have been confirmed 
by the available test material. 

From this comparison of the various experiments, we obtained 
1 two important indications for interpreting oil film photographs 
and Schlieren photographs. In oil film photographs, the vortex 
center has to be assumed above the separation line, so that the 
typical oil lines only show half of the vortex flow. This is in 
contrast to the evaluation method of Rao and Whitehead 15^, which 
places the vortex core exactly in the center between the separation 
line and the attached flow. In Schlieren photographs, we can ob- 
serve a similar condition. There the dark gray and not-sharply- 
defined regions only are indications of the upper part of the vor- 
tex to the core. The lower vortex region in the vicinity of the 
wing surface can no longer be seen because of interferences and the 
fact that the body surface is so close. These observations about /29 
the vortex position in oil film photographs and Schlieren photo- 
graphs are also supported by investigations of Mannerie and Werle 
175J, as well as of Drougge and Larson Il4]. 
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Prom the result for M^ « 2.5 and a * 10°, we can now investi- 
gate the Influence of angle of attack for the same incident Mach 
number on the flow field. We can expect that leading edge separa- 
tion in the sense defined above will occur, and only the geometry 
of the flow field will change. The Schlieren observations for 
angles of attack between a = 10° and 25° allow one to analyze 
the relationship with the Pitot isobars over the wing cross-section. 
Figure 28 shows that up to a = 15°, the primary vortex separates 
from the surface, whereas the secondary vortex essentially keeps 
the same position. The region Influenced by the vortices expands 
further into the direction of the symmetry plane and there reduces 
the region of attached flow. This development of the flow field 
up to average angles of attack is normal in the regions of leading 
edge separation and was to be expected. Prom Pig. 20b we now find 
that for higher angles of attack after a, = 17. 5 ° » there are some- 
times substantial deviations from the initial flow conditions. 

The primary vortices move away further from the surface and then 
increase their area of influence, so that no more attached flow 
occurs in the central region. The vortices on both wing halves 
are close to one another in the symmetry plane and in this way in- 
duce a very large velocity component in the direction of the top 
side. The flow, however, must again be parallel to the surface in 
the vicinity of the lee side and therefore is deflected through a 
shock which can only be found in the direct surroundings of the 
central plane. As can be seen from the Pitot isobars, for example, 
at a = 20°, this central plane shock is approximately parallel to 
the surface of the wing. When the shock occurs or when the vortices 
meet in the symmetry plane, the secondary vortex could not be found 
either in oil film photographs and Schlieren photographs or in the 
Pitot pressure distributions. A dead-water region is formed in the 
outer wing area which is shown as a shaded region in the Pitot iso- 
bars . 
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The continuity of the flow for transition from small angles 
of attack to large angles of attack is confirmed by the pressure 
distributions over the span in Fig. 21. From the curves, we show 
the positions determined from oil film photographs for attached 
flow in the central region (A), secondary vortex (SZ) and primary 
vortex center (PZ). Up to a = 10°, one observes typical pressure 
variations for leading edge separation with a constant central 
region and a strong pressure drop in the region of the vortex 
center. This result was also obtained from the CUED wind tunnel 
and shows good agreement with the ILR tests considering the dif- 
ferent tunnels and measurement methods. For the higher angle of 
attack range, these clear statements can no longer be made. The 
entire pressure level lies very close to the vacuum limit, so that 
the effect of vortex flows on the pressure distribution is limited 
by this. (The pressure increase in the outer wing area at a = 25° 
is already an influence of the disturbance which comes from the 
tall region, recognized in Schlieren photographs.) The result 
from oil film photographs can be compared very well with the pres- 
sure distribution and support previous results and models . 

Figure 22 shows the result of oil film photographs in Schlieren 

photographs and gives the positions of the attached flow in the 

central region and the position of the vortices above the wing for 

= 2 . Q and 2.5. The typical change in the position of leading 

edge vortices with increasing angle of attack is known from tests 

by Monnerie and Werle 75 and the theoretical analysis of Pershing 

32 . The agreement with the present vortex positions determined 

and these results show’ that for M = 2.0 we can draw the conclu- 

00 

sion that there is leading edge separation for the angle of attack 
region ‘ up to a = 20°. For M ro = 2.5, on the other hand, above 
a = 15° there is an unexpected change in the vortex position, which 
indicates that then there is no classical leading edge 

separation. 


Summarizing, we can say that to the left of the Stanbrook- 
Squlre region between a ■ 5 0 and 15° there Is leading edge separa- 
tion with primary and secondary vortices. Prom a further Increase 
of the angle of attack, a flow field results which differs very 
greatly from the previous Ideas, which show a dead-water region 
In the region of the former secondary vortex and which has a shock 
In the symmetry plane parallel to the surface. Since similar 
changes also occur for M^ = 2.0 but at higher angles of attack /31 

(a > 20°), here we have to specify a new limit. This result Is 
given here In this paper only as an Indication, because the pur- 
pose of our Investigation Is to give a detailed analysis of the 
flow regions Inside and to the right of the Stanbrook-Sqiiire bound- 
ary. 


4.2.2. The Stanbrook-Squlre region 

Since the flow regions to the left of the Stanbrook-Squlre 
region have now been discussed and the typical characteristics 
of leading edge separation have been mentioned, we will now in- 
vestigate the transition to shock-induced separation. For this 
purpose we will again assume the same incident flow conditions 
with the angle of attack of a = 10° and the Mach number 1^, = 2.5- 

The pressure coefficients on the flat lee side are shown in 
Pig. 23 over the span for different Mach numbers but for a fixed 
angle of attack a = 10°. One can see large differences in the 
pressure distributions only along the outer wing region where the 
vortex dominates. The pronounced underpressure peak, already dis- 
cussed somewhat for M^ = 2.5, becomes flatter for higher Mach num- 
bers until finally the entire region only has a constant pressure. 
Even at higher Mach numbers outside of the St anb rook-Squire limit, 
the form of the pressure distribution no longer changes but the 
pressure level only is displaced to higher values. In the central 
region where attached flow is assumed, the position of the separa- 
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tlon point (A) only changes slightly. The pressure In this region 
is constant and Is only changed slightly with a change in Mach 
number. An additional evaluation of the pressure distribution 
is given in Pig. 24, Here the pressure In the region of attached 
flow, as well as the minimum pressure coefficient (vortex center), 
is shown for Increasing normal Mach numbers M N - In the region of 
the St anb rook-Squlre boundary, the pressure coefficients reach a 
maximum in the central plane. The pressure in the vicinity of the 
vortex center decreases ahead of and behind the Stanbrook-Squire 
boundary more steeply than inside of this region. In addition, we 
show the corresponding values for the delta-shaped lee side, which 
can be correlated with the other results, at least for the pressure 
near the symmetry plane. 

Prom the pressure distribution we can see that the central 
region with the parallel flow remains unchanged with increasing 
Mach number, which is indicated by the combined oil film photo- 
graphs and Schlieren photographs of Pig. 25 . The vortex center, on 
the other hand, is displaced towards the central plane so that the 
diameter of the vortex must become correspondingly smaller. In 
addition, these figures show that inside the Stanbrook-Squire bound- 
ary, the Intensity of the secondary vortex becomes smaller and it 
is hardly perceptible at M ro = 2.8. When the secondary vortex 
vanishes, only a dead-water region exists between the leading edge 
and the primary vortex. 

In addition, the Schlieren photographs show the approach of 
the leading edge shock to the wing. In this way the Mach number 
component Mjj reaches values above one. It is to be expected that 
immediately downstream of the leading edge there will be supersonic 
conditions above the dead-water region and the character of the 
lee side flow is changed in this way. 

Two Schlieren photographs of the flow at a = 10° and different 
Mach numbers are given in Pig. 2 6. The dark regions which can be 
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Interpreted as the upper halves of the vortices on the lee side 
are quite sharply limited first at M^ » 3.0 in the direction of 
the wing surface, but become more diffused at higher speeds. 

Two dark lines emerge from the trailing edge at the lee side of 
the model, which have to be the result of different wing regions. 
The upper region Is in the form of a fan and can be Interpreted 
as an expansion of the flow in the central region at the trailing 
edge. The second narrow line is almost parallel to the top side 
and can be interpreted as a wake dead water from the leading edge 
region of the wing. 

Figure 27 shows an evaluation of the Schlleren and oil film 
photographs as a function of Mach number and gives the position 
of the vortices. We should mention that up to the present only 
little information about the Mach number dependence of the vortex 
position is available. The only systematic investigation of 
Okerbloom and Sarantsev [16] showed an almost constant height of 
the vortex above the wing between M ro = 2.0 and 4.0. This result 
was not confirmed here but one can observe that the vortex moves 
closer to the top side, which agrees with Monnerie and Werle [ 75J • 
Only when the Stanbrook-Squire region is crossed does the position 
become fixed for a fixed angle of attack. 

Summarizing, we can establish the following for the transi- 
tion region (Stanbrook-Squire region). The present results for 
an angle of attack a = 10° and different Mach numbers confirm a 
change in the type of flow above the Stanbrook-Squire region. 

This is a continuous process which is characterized by the follow- 
ing flow behavior: At low incident Mach numbers there is clearly 

a leading edge separation with primary and secondary vortices, as 
well as a parallel flow in the region of the symmetry plane. The 
increase in the incident speed and the increase in the normal com- 
ponent with respect to the leading edge means that the higher 
energy flow can resist the separation at the leading edge longer 


and then a deflection around the leading edge becomes possible. 

In this way the separation angle between the surface and the 
feeding vortex surface becomes smaller and the entire vortex be- 
comes flatter. In this way the diameter of the primary vortex 
becomes smaller and therefore the region of attached flow near 
the central plane can remain constant for the most part. Because 
the height as well as the Intensity of the primary vortex de- 
creases with Increasing Mach number, the velocity components 
which trigger the secondary vortices are greatly reduced. About 
in the center of the S-S region (M^ » 3*0), the secondary vortex 
can no longer exist and the region between the leading edge and 
the primary vortex has to be Interpreted as a dead-water region. 

In addition to the flow conditions in these regions with high 
friction, one can also observe changes in the external flow. In 
the investigated models, there is a separated leading edge shock, 
but within the S-S region the speed of sound Is reached at the 
leading edge. Above the layers with friction, a supersonic ex- 
pansion occurs In the vicinity of the leading edge which determines 
the pressure level in this region. 

The greater the Incident flow speed, the flatter the layer /34 

with friction becomes over which the expansion occurs. In this 
way the expansion is intensified and there is a substantial deflec- 
tion of the flow in the direction of the central plane. The symmetry 
condition in the central plane can be satisfied by the deflected flow 
only by means of a shock. During the increase in Mach number, the 
shock intensity increases and the vortex intensity decreases. 

Summarizing, we can say that the following events characterize 
the transition region for constant angle of attack and increasing 
Mach number: 

— The vortex is flatter above the surface and reduces its 
influence region whereas at the same time the intensity 
of the secondary vortices is decreased. 
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— When the secondary vortex no longer occurs, a dead- j 

water region Is formed Inside the layer with friction. j 

j 

— The flow in the leading edge region reaches supersonic 
conditions and the starting expansion over the friction 
layer results in an Internal shock. 

— Even though a continuous transition occurs and therefore 
one cannot give a fixed boundary for the change in the 
flow types, at least for M * 3.0 we established that 
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this results In the most noticeable changes in experl- I 

j 

ments . j 

A detailed analysis of the flow conditions to the right of 
the St anb rook-Squire region has to be added to the above and will 
therefore be carried out In the following section. 

4.2.3* Separation with a shock 

In the following we give an example of an analysis of the 
flow field in the region of the so-called shock-induced separa- 
tion for a = 10° and M ot = 3.5 and in a discussion of the present 
experiments. Figure 28 shows different test results for the 
multiple comparisons. From this we can see that a dead-water 
region exists directly on the surface in the outer wing region. 

Between it and the external flow there Is a type of boundary 
layer which is probably highly turbulent. In the central region, 
there Is attached flow which is aligned parallel to the symmetry 
plane . 

j Figure 29 gives an additional evaluation of Pitot pressure /35 

measurements and also an additional model. We show the Pitot 
pressure ratios over height for different span positions. In 
these experiments one can assume that the Pitot probes aligned 
parallel to the incident flow will provide reliable measurements 
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because of the high Incident flow speed and the slender models, 
because In this case the transverse speeds are of secondary Im- 
portance. This means that Pitot pressure ratios of V p p 0o *' 
represent conditions like in the region of undisturbed Incident 
flow or If the value Is zero, they can be Interpreted as dead- 
water regions. Figure 29 gives an evaluation in the flow model. 

Based on the experimental results, we can establish the 
following model. At the leading edge there Is flow separation 
just like before, but its intensity and momentum is very small 
and therefore can be called dead water. Since the separation is 
assumed to be a turbulent boundary layer, one can count on a re- 
attachment (according to the separation bubble in profile theory) . 

On the other hand, the processes in the outer flow also influence 
the boundary layer. The supersonic flow which exists at the lead- 
ing edge expands over the contour of the separated layers and for 
symmetry reasons, has to be deflected backwards again by means of 
a shock. The pressure jump through the inner shock affects the 
boundary layer and results in a reattachment of the separated flow. 
Because of this forced premature reattachment, according to profile 
theory we can speak of a separation bubble. Therefore, a stagna- 
tion results, which separates the dead-water region from the 
attached flow in the central region. Within the separation bubble 
there is a reverse flow which occurs, which has a certain circula- 
tion and which again can be called a vortex. 

Therefore, we have established that in contrast to the original 
assumption, there is no shock-induced separation. In the case of 
the wings investigated here, the flow separation is not produced by 
a compression shock, but is a leading edge separation and the sepa- 
rated flow is made to reattach because of the shock in the external 
flow. It is quite possible, as found for small angles of attack, 
that a separation can occur at the leading edge without a shock. /36 

On the other hand, just like for high Mach numbers, a supersonic 
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expansion at the leading edge and a deflection over a weak Internal 
shock is thinkable without having a separation occur. The wings 
which are being discussed have been investigated in exactly this 
state range, where there is both separation at the leading edge as 
well as supersonic expansion with subsequent Internal shock. Both 
flow components occur Independent of one another, but Influence 
each other until an equilibrium state is established between them. 

Thei’efore, we have shown that for the wings under discussion, 
the Stanbrook-Squire region is only a boundary with respect to the 
flow in ths frictionless region, because to the left and the right 
of the boundary there is a separation at the leading edge. The 
region to the right of the transition region will be called "separa- 
tion with shock" region according to the results discussed. 

Based on the previously developed flow model, we will now in- 
vestigate the changes which occur for a variation of the angle of 
attack if we assume a fixed incident Mach number = 3-5. As 
shown from Pig. 2, experiments are located in the region which up 
to the present was called "shock-induced separation" region. The 
oil film photographs in Pig. 30 show the development of the lee 
side flow for = 3*5 fo’” small-to-medium-sized angles of attack. 
Por angles of a = 1° the flow field appears to be non-conical and 
is quite similar to the experiments of Rein [46] , and its flow 
model shown in Fig. 6. A second interpretation of these oil film 
photographs is that the transition from laminar to turbulent bound- 
ary layer causes the observed differences. However, for each of 
these possibilities, we can assume a vortex formation at the lead- 
ing edge, until an angle of attack of a = 5.5° is reached. Por 
this angle of attack, the oil film photographs already have charac- 
teristics which were found when crossing the S-S boundary. Por 
a - 8°, we again have the typical situation of a separation bubble- 
with an internal shock, which exists up to angles of attack of 
a = 16.5°. However, at these angles there is already a great in- 
fluence of the tail region, so that only the region near the tip 
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has to be considered as not being Influenced. 


The observations of the oil film photographs made are well- 
confirmed by the pressure distributions In Fig. 31. The vortex 
positions (WZ) taken from the oil film photographs and the regions 
of attached flow In the central plane (A) are shown on the corre- 
sponding curves. For small angles (a * 1° - 3.5°), one obtains 
pressure variations along the leading edge which are typical for 
vortex formation. At a = 5.5° the transition to the solution 
bubble with internal shock seems to appear. The discrepancy In 
the pressure distribution and the oil film photograph can be ex- 
plained by the testing technique, because In all experiments the 
models have almost zero angle of attack when the wind tunnel Is 
started and later on the desired angle Is regulated using the 
pneumatic system, l.e.. In the start-up position certain oil lines 
already form which will hardly change afterwards . At higher angles 
of attack, again the influence of the tail becomes noticeable, 
which at a = 12° produces a pressure increase first In the outer 
wing area. This disturbance from the wake Is Investigated in a 
special chapter as far as it is relevant for the present experi- 
ments, and will be critically analyzed. 

Figure 32 shows the vortex positions for different angles of 
attack. The position of the attached flow In the central region 
is also shown and only has a large change for small angles. The 
further development decreases continuously up to the angle of 
attack a = 13°, where the vortex extends to the symmetry plane. 

The position of the vortex is displaced between a = 3. 5° (or 
5.5°) and 8° very greatly towards the central plane, and when 
there is a further increase in the angle of attack, its position 
remains almost unchanged. This Indicates a different flow field. 
For small angles of attack, the position in the span direction is 
comparable with incompressible (M m =0) results, so that here vor- 
tex formation at the leading edge can be assumed. This was con- 




firmed by experiments of Squire [4l] , as well as Monnerle and 
Werle 1753, In which both the pressure distribution as well as 
the oil film photographs show the same characteristics, as occur 
here when the angle of attack is Increased between a = 4° and 8°. 

For a = 8° we can already expect the separation bubble with 
the internal shock. Up to a = 12° there is only parallel attached 
flow In the central plane. At higher angles of attack, the vortex 
center again is displaced in the direction of the leading edge and 
also separates from the surface, without the occurrence of parallel 
flow in the symmetry plane. The vortex therefore expands its /38 

region of influence and apparently at a = 20° again takes on values 
according to an approximation from the available data as is given 
by incompressible flow. This tendency of a separation line moving 
downwards was also demonstrated by Cross [491, who established 
leading edge separation for angles of attack of a > 20° for rela- 
tively thin delta wings at M ro = 10. However, this contrast with 
the static pressure measurements in which the pressure increase 
measured near the central plane is displaced with increasing angle 
of attack in the direction of the symmetry plane. On the other 
hand, it is quite possible that the internal shock in the external 
flow moves somewhat inwards and the related pressure is felt up to 
the top side. At the same time, in this angle of attack region 
the close vacuum again influences the pressure level. 

We already described how one can derive information about the 
vortex intensity from the oil film photographs . Figure 33 shows 
the relative flow angle for the investigated Mach number range 
which characterizes the vortex intensity, and it is plotted as a 
function of angle of attack. For all of the regions to the right 
and left of the S-S region, we can observe a maximum in intensity 
which is related to a meeting of the vortices in the symmetry plane 
and the related displacement of the attached flow. For Mach num- 
bers up to the center of the S-S boundary (M^, ~ 3*0), these intensity 
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maxima are equivalent to the last occurrence of the shock in the 
symmetry plane. This central line shock at higher Mach numbers 
can no longer occur, and it is caused by the change in the flow 
type. This is because the intensity of the vortex and therefore 
the velocity component towards the surface is too small. A com- 
parison of the vortex Intensities for the different Incident Mach 
numbers allows one to guess that at about M = 2 the maximum vor- 
tex intensity is reached. 

Summarizing, we can describe the flow type "separation with 
shock" as follows for the case of a varying angle of attack: 

— For small angles of attack the vortex separation from 
the leading edge could be observed (a = 0° to 3.5°). 

— After a transition region similar to the S-S region, 
there is a separation bubble at the leading edge and 
an internal shock (a = 8° to 14°) up to medium-sized 
angles of attack. 

— For high angles of attack we found an enlargement of 
the vortex and a new vortex formation from the leading 
edge. However, clear data are not available in this 
case, because the strong tail influences falsify the 
investigations (a < l 1 ! 0 ). 

4.2.4. Classification of the flow types 

In the previous sections, we investigated the flow region for 
the flat leeward side, which according to Fig. 17 includes the 
Stanbrook-Squire region. The consequences drawn from this and the 
models developed are shown in Fig. 34 in a - M N diagram. In 
contrast to the previous S-S region, the region differs for small 
angles of attack. Here there is leading edge separation also on 


the right side of the S-S region, and the new transition region 
approaches the M N axis with Increasing Mach number. Accordingly, 
there Is also an upper boundary at which the transition from the 
separation bubble with Internal shock to the leading edge vortex 
formation occurs. However, here we can also assume that above 
the vortices, shocks are produced. These models can be derived 
from the extreme case a = 90°, where the dead water is closed off 
by shocks, similar to a tail flow. 

Comparisons between the discontinuities in the "thin shock 
layer" theory and experiments made Squire [6l ] assume that for 
thicker wings, the boundary between the regions B and D is dis- 
placed to smaller angles of attack a N . With the present results, 
we can draw the conclusion that with increasing wing thickness, 
i.e., greater angle between the top side and the bottom side per- 
pendicular to the leading edge, the St anb rook-Squire boundary 
runs into the region B for large and small angles of attack, as 
is also shown in Pig. 3^ . The region B of the "shock-induced 
vortex formation" is limited on the left side for all wing thick- 
nesses primarily by the Mach number, but for thick wings it is 
also limited upwards and downwards by the angle of attack. 

Summarizing, we can therefore delimit the complex lee side 
flow field, not only by the Stanbrook-Squire region but also by 
regions which depend primarily on a or M a . 

4.3. The influence of the cross-section shape on the flow 
for separations with internal shock 

In the previous discussion, the flow field of the delta wing 
was treated with a flat (,'plane) leeward side. The bottom side of 
the wing was delta- shaped. In the following, we will consider 
delta wings with non-flat (delta shaped or conical) top side. The 
purpose of this investigation is to obtain Information about the 
influence of the cross-section shape on the lee side flow. The ex- 


perlments are restricted to the range "separation bubble with in- 
ternal shock," i.e., to the region to the right of the Stanbrook- 
Squlre region. In most results It was assumed that the bottom 
and top side can be treated separately In the analysis. In a 
further section we will briefly discuss the influence of the 
bottom side on the lee side flow. 


4.3.1. The delta shaped lee side 

The delta shaped lee side was selected as the representative 
top side for investigating the influence of geometry. The shape 
was used primarily in the few known experiments, so that there is 
some capability of making comparisons. For a delta-shaped lee 
side, there are two plane model halves and therefore it can be 
assumed that only the angle 0 Cbetween the leading edge plane and 
the ridge line) is a geometric parameter which has any influence. 

In the description of the flow fields for the delta-shaped 
lee side, we select an incident state with M^ = 2.5 and a = 10° as 
the initial example. After this we will demonstrate the influence 
by changing the angle of attack and increasing Mach number. 

Figure 35 shows a flow model for the initial case which is the 
result of experiments. In the plane perpendicular to the leading 
edge, we find an angle between the incident flow and the top side 
of oijj = 17°. The angle of attack is therefore about half as 
large as was was the case for the corresponding incident flow 
states with the flat lee side. Among other things, this means that 
at the leading edge there is a Prandtl-Meyer expansion with attached 
flow on the top side, i.e., in the region of the leading edge there 
is no separation. The flow is deflected by means of the Prandtl- 
Meyer expansion, however, so that the resulting flow direction goes 
towards the symmetry plane. An internal shock, i.e., a series of /4l 
compression waves, provides the reverse deflection which is required 
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to satisfy the symmetry condition. The compression, however, re- 
sults in a separation with the vortex. These processes are re- 
stricted to the outer wing regions and the vortex has a relatively 
small size, so that again there is attached flow In the vicinity 
of the ridge line. 

This model is only partly applicable for the flow, if the 
angle of attack is changed together with a constant Mach number. 
Figure 36 shows the results of Pitot pressure measurements for 
Moo ■ 2.5 and different angles of attack, and” it also shows the 1 
results of oil film photographs and Schlleren photographs. For 
small angles of attack up to a = 7°, there is leading edge vortex 
formation. Their feeding vortex surfaces approach the surfaces 
with increasing angle of attack and finally they attach. Between 
a = 7° and 11°, the normal component of the incident speed is so 
large that there is no leading edge separation and the attached 
flow exists in the outer wing region. Up to an angle of attack 
of a = l4°, we can therefore assume a similar flow, as already de- 
scribed for a = 10°. When the angle of attack is increased further, 
the vortex center on the one hand moves in the direction of the 
central plane and, on the other hand, moves away from the wing sur- 
face. Since the internal shock does not change its position so 
much, a dead-water region forms which increases inside between the 
separation line and the vortex center. The vortices of the two 
wing halves are then so close to the ridge line, that there is no 
attached flow in the central region anymore. 

The Mach number influence on the lee side flow for a delta- 
shaped surface and fixed angle of attack a = 1Q° will be discussed 
using the pressure distribution shown in Fig. 37- An increase in 
the Mach number here results in a correspondingly larger normal 
component M N , which is decisive for the expansion at the leading 
edge. Based on the greater expansion rate, the internal shock and 
therefore the separation caused by it are displaced in the direction 
of the central plane. These processes can be observed for Mach 


numbers of up to M ro = 3.0. 


The flow conditions which prevail at higher Mach numbers 
can best be explained using the model shown In Fig. 38. The 
most characteristic feature Is the oil film photograph shown 
there, which shows that the flow is only partially conical. 

Here there Is a certain similarity with the flow around the flat 
lee side at M m = 3*5 in the lower angle of attack range. There, 
there was a partial non-conlcal flow field (Fig. 30) along the 
flat leeward side in the tip region of the wing, apparently 
caused by the transition of the boundary layer. This transition 
from laminar to turbulent boundary layer Is also Influencing the 
states on the top side in the case of the delta-shaped lee side. 

In the conical tip region, just like before, one finds an expan- 
sion in the region of the leading edge with a subsequent internal 
shock, which brings about a separation of the laminar boundary 
layer. After about 30 % of the wing chord, the transition to the 
turbulent boundary layer starts. At the leading edge, a short 
separation bubble is formed. By means of this change, the vortex 
which occurs in the tip region of the wing loses its ’'supplying" 
vortex surface and exists as a free vortex downstream of the transi- 
tion region. This region of the wing has to be considered non- 
conical. 

For a critical evaluation of the flow model, we should indi- 
cate that for =3.5, we carried out the greatest number of tests 
with bodies having 0 =^= 0° in the CUED wind tunnel. It is known 
that the transition from the laminar to turbulent boundary layer 
is influenced substantially by the Reynolds number, which in the 
CUED wind tunnel is about 5 times greater in the ILR tunnel. This 
aspect has to be considered for the flow field analysis carried 
out here, especially for transition to the higher Mach numbers. 
Therefore, in the following discussion we will separately discuss 
the tests in each wind tunnel, in order not to have to consider 


Re number effects which also occur In this flow range. 


Starting with the results for M w * 3.5 and a ■ 10°, using 
the pressure distributions of Pig. 39, we Investigate the in- * 
fluence of the angle of attack for a fixed Mach number. The out- 
line pressure distributions are complemented by information ob- /*>3 

talned from oil film photographs about the position of the 
attached flow (A), the vortex center (PZ) and the separation 
bubble (B). At a * 0° the flow is attached along the entire top 
side and only near the leading edge does a weak expansion influence 
the pressure values. Already at a * 3° one can notice a conical 
vortex pair in the tip region of the wing, whose intensity, however, 
is so weak that after transition in the central region there are 
almost no more influences anymore. With increasing angle of attack 
the intensity of the free vortex increases, but at the same time 
its region of influence decreases. The separation bubble (B) at 
the leading edge, on the other hand, becomes greater for higher 
angles of attack, whereas the internal shock and therefore the 
region of the attached flow is displaced further towards the central 
plane . 

The span positions of the vortex center and the attached flow 
are shown in Pig. 40 for all angles of attack and Mach numbers. 

Prom the vortex positions and by means of a comparison with the 
values for the flat lee side (JV^, = Q and 2.5), we can see that for 
the delta-shaped lee side for small Mach numbers and small angles 
of attack, there is also vortex formation at the leading edge. The 
position of the vortices develops quite differently as soon as 
supersonic expansion occurs at the leading edge and therefore the 
internal shock and separations occur with vortex formation. It 
seems possible to find a limit which characterizes the transition 
between the flow types,. This region is shown in the diagram by a 
shaded surface and shows good agreement with the transition regions 
shown in Pig. 3^ from the flow around flat delta wings with small 
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angles of attack. The position of the vortices Is relatively In- 
dependent of the He number Influence because it Is always created 
in the vicinity of the model tip. When we consider the attached 
flow, we again have to consider the different wind tunnels. 

Summarizing, we can say the following about the flow field 
of the delta-shaped lee side In the investigated Mach number range: 

— In contrast to the flat top side and caused by the shape 
of the lee side, small angles of attack and large Mach 
numbers M N occur perpendicular to the leading edge. Prom 
this, attached flow results at the leading edge, together 
with a supersonic expansion. The Internal shock Induces 
a separation. 

— For small angles of attack, there is again vortex formation 
at the leading edge. 

— For certain incident flow conditions CM -, Re-influence), 
there can be a change in the flow field described above. 
Near the wing tip there is a transition from laminar to 
turhulent boundary layer. Downstream of this region, a 
short separation bubble is created which is closed off by 
the internal shock. The separation produced in the tip 
area exists as a free vortex above the wing. 


4.3.2. The influence of different shape parameters (angle 6) 

In order to describe the influence of various top side shapes 
on the flow field of the lee side, we first assume that the differ- 
ences in the cross-section shapes are essentially described by the 
angle 0 (between the leading edge plane and the ridge line). 

Four different lee side shapes with angles between 0=0° and 


14 U were available. This means we can describe the transition of 
the flow fields from the flat lee side to the delta shaped lee 
side. Figure 41 shows oil film photographs of the top sides for 

M • 3*5 and a » 10°, and the corresponding pressure distributions 
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are also shown. With increasing angle ©, the normal Mach number 
M n becomes larger, whereas at the same time the normal angle of 
attack a N becomes smaller. From the pressure distributions, we 
can see that the pressure level also decreases with increasing 6. 

The values shown from the Prandtl-Meyer expansion with the cross - 
flow components around the leading edge confirm this tendency. 
Quantitatively, however, there are substantial differences between 
these calculated pressure values and the measured pressure values, 
which among other things are due to the influence of the leading 
edge shock. All of the investigated wings have differently shaped 
lower sides and therefore give a different position and shape to 
the shock. In this way, behind the shock ^different values for 
M^j and are created than calculated for the incident flow in 
the normal plane. 

Figure 42 shows the influence of the shape angle 0 on the flow 
field for different top side shapes. Both the positions of the 
vortex centers as well as the attachment line vary slightly linearly 
with increasing 0 angle in the direction of the central plane. We 
should note that a reduced. semi-span has been introduced which is 
the ratio of the lee side span for flat (0 = 0°) and super- 
elevated wing (e > 0°). 
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For the transition from a flat to a delta-shaped lee side, we 
can develop the following model for the flow processes. Starting 
with the flow model with separation bubble at the leading edge. 


Internal shock and attached flow In the central region (Pig. 39), 
the Mach number In the normal plane becomes greater as the 
shape angle 0 Increases, whereas the angle decreases. The 
flow therefore has more energy and the separation bubble at the 
leading edge becomes flatter. In this way the Intensity of the 
circulation Inside the separation Is changed, so that the radial 
velocities decrease compared with the velocities in the axial 
direction. The reattachment line of the separation bubble Is 
displaced toward the central plane and the region of the attached 
flow Is reduced In size. Using the previous flow models, we must 
conclude that the Internal shock changes its position only slight- 
ly. For sufficiently high transverse speeds, the transition of 
the boundary layer starts and determines the processes on the lee 
side already discussed In the previous chapter. The boundary layer 
transition is such that already for 0=9° it can be observed along 
the rear part of the wing. 

The following can be stated about the influence of the shape 
parameter 0 for plane top sides : 

— With Increasing angle 0 the separation bubble becomes 
flatter and moves in the direction of the central plane. 

The Intensity of vortex motion decreases. The inner 
shock harely changes its position. 

— For higher angle © Incident flow states are reached 

which can bring about a transition of the boundary layer. 

— The geometric changes of the flow field are linear In the 
shape parameter 0. 

— For the same incident flow conditions and M^, we observe 
basic changes in the flow field for the top sides, with 
different angle 0 (within the boundaries mentioned here). 
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The present experiments were all carried out In the range of 
flow with the separation bubble and Internal shock. For the In- 
vestigation of these wings Inside the Stanbrook-Squlre region, we 
can therefore no longer assume a linear dependence on the shape 
parameter. The boundaries In the <* N - diagram therefore can only 
be assumed to be valid for a single angle 

4.3.3* The Influence of different leading edge shapes (angle ¥ ) 

In the previous section we discussed the effects of different 
top side shapes on the lee side flow field In such a manner that 
we found a direct dependence on the shape angle e. In the follow- 
ing we will fix the angle 0 and discuss the Influence of different 
leading edge angles ¥ (measured perpendicular to the leading edge 
between the top side and the bottom side), l.e., different curva- 
tures of the lee side. Two wind tunnel models were used for this: 
the cone model (Fig. 43) and the modified cone model (Fig. 44). 

The results are compared with those for the delta model because all 
three wings have the same shape angle ('0 * 14°). The leading edges 
of the Individual wings are shaped In such a manner that the conical 
delta wing has the largest angle (¥ = 75°) and the modified cone 
model has the smallest angle (¥ = 18°), whereas In the delta model 
we measured ¥ = 40°. 

The delta wings described here all have flat undersides and 
the plane of the leading edge Is the same as the lower side plane. 

If one assumes a delta model with a plane lee side half, then the 
angle of attack in the plane perpendicular to the leading edge Is 
described by the angle between the normal speed and the lee side. 
Therefore, for « = 0° we also have = ¥ . If we have curved sur- 
faces for a fixed shape angle 0, then we also find different angles 
between the normal speed and the tangential plane in the region of 
the leading edge. 


50 


Very large leading edge angles f, as are found for the cone 
model, result; In a flow In the region of the leading edge where 
the lee side Is seen under a negative angle. In this way, as /47 

shown in Pig. 43, a stagnation point Is formed on the top side. 
Downstream of It, the flow field found from Pitot Isobars and 
oil film photographs show mostly agreement with the flows along 
the lee side of circular cones with angles of attack. Prom these 
investigations, we can see that at least for small and medium 
angles of attack the leading edge no longer belongs to the lee 
side. In contrast to the delta wings used up to the present with 
a sharp leading edge and therefore a fixed separation line, the 
leading edge is of secondary importance for the lee side. 

The other extreme case is a very small leading edge angle ¥, 
as for example occurs for the modified cone model (Pig. 12). In 
the vicinity of the leading edge, the cross-section shape Is coni- 
cal with the angle Y » 18°, whereas the central part consists of 
a cone segment. Here again the size and direction of the Incident 
flow component perpendicular to the leading edge Is unchanged com- 
pared with the delta shaped lee side, so that for an angle of attack 
of a = 10° an expansion angle of a Nex * 38° is achieved in the lead- 
ing edge region. As an example, Pig. 44 shows the pressure distribu- 
tion and the characteristic flow model for two angles of attack. 

For angles of attack around a = 10°, there Is an expansion at the 
leading edge which induces a low pressure. Downstream the flow is 
uniform and the pressure distribution is constant. Only a very 
strong change In the flow direction caused by the body contour re- 
sults in an Internal shock whose Intensity is greater than that 
found for previously observed delta wings. In the conical central 
region, we can again observe the vortex system which is typical for 
this shape. For angles of attack of a > 15°, we can only achieve 
an incomplete deflection and expansion of the flow, so that total 
separation starting at the leading edges occurs. Inside of the 
separated region, there is subsonic flow and neither the internal 
shock nor other Influences change the distribution here, which is 
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approximately constant over the wing. For higher Mach numbers, 
this flow field can already be observed for smaller angles of 
attack than a = 15°. 

From these examples we already can see that for a thick delta 
wing with a fixed shape parameter 0 the angle ¥ between the top 
side and the bottom side perpendicular to the leading edge sub- 
stantially influences the aerodynamic properties such as lift, drag 
and efficiency of the control elements. Excessively large angles 
can lead to stagnation points in the leading edge region and there- 
fore to overpressures on the lee side. On the other hand, the flow /^8 
for small angles ¥ and therefore large flow angles of the leading 
edge have a tendency to incomplete expansion or to complete separa- 
tion at the leading edge. 


4.3.4. The influence of the cross-section shape of the lower side 

The lee side flow can only be influenced by the underside if 
the leading edge shock is separated. This applies for all delta 
wings in the investigated measurement range and means that it is 
necessary to estimate the influence of the lower side shape on the 
lee side flow. The position of the leading edge shock is primarily 
determined by the shape of the underside. Therefore, the inclination 
of the shock near the leading edge is, for the most part, determined 
by the shock distance in the vicinity of the central plane of the 
delta wing. Since the separated shock causes a deflection of the 
incident flow, different inclination angles of the leading edge 
shock can lead to substantial changes in the flow field on the lee 
side. 

Figure 45 shows the different shock shapes of concave (wave 
rider) and flat (delta model) undersides for the identical lee side 
shape. At the same time, the Pitot isobars allow an analysis of 
the flow field on the top side for different shapes of the leading 
edge shock. For small angles of attack (a w 9°) and for the flat 
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delta model, there Is already a separation with an Internal shock 
and vortex formation above the wing. On the other hand, the wave 
rider has a completely different behavior. On Its top side, there 
is a type of separation bubble In the vicinity of the leading edge, 
and there Is attached flow over large parts of the wing. This 
discrepancy must result primarily from the position and shape of 
the leading edge shock. For a flat underside, this lies both in 
the plane of the leading edge as well as In t' ’• symmetry plane at 
the same distance from the wing. The shock is almost perpendicular 
to the leading edge plane in the region of the leading edge plane 
and therefore brings about a deflection of the flow. The deflec- 
tion angle at the leading edge becomes larger, whereas the corre- 
sponding local Mach number is reduced but remains in the supersonic 
range. For the wave rider, the shock is relatively close to the 
plane of the leading edges and is barely curved. The flow compo- 
nent perpendicular to the leading edge therefore reaches an almost 
perpendicular shock and retains its initial direction, approximate- 
ly. At the leading edge, there is no expansion but a separation 
bubble is formed. For average angles of attack (a = 1*1°), the 
shock for the wave rider is substantially more curved in the region 
of the leading edges, so that we find a similar flow picture as 
was found for the lee side of the delta model. Nevertheless, be- 
cause of the separation of the shock from the leading edges, we 
still find differences which are characterized by the vortex posi- 
tion and the size of the attached flow in the central region. 

Theoretical results of Squire [6l] represent a further indica- 
tion about the relationship between the flow fields on the top side 
and the bottom side. The calculations showed that on the underside 
of thin delta wings, there is a sudden transition from a stagnation 
line in the symmetry plane to a second one which is near the lead- 
ing edge. This change only occurs for certain incidence angle con- 
ditions and from comparisons with known experiments we can correlate 
this result with the transition from leading edge separation to 


shock-induced separation on the lee side. These theoretical re- 
sults are supported in the following by means of an example for 
the delta model. Assuming that the change In the flow type on the 
lee side is directly related to the change from the case with a 
stagnation line In the symmetry plane to two stagnation lines on 
the lower side, the change has to be continuous In contrast to 
the theoretical results. This assumption is supported by Pig. 46, 
which gives the stagnation line positions on the flat pressure 
side as a function of Mach number. According to this, because of 
symmetry reasons, a stagnation line in the central plane is ob- 
tained and, depending on the incident flow conditions, there is an 
additional one in the region of the leading edge. For the region 
of separation with internal shock, we have two stagnation lines. 

For attached leading edge shock as well as in the region of the 
leading edge vox’tex, there is only the symmetry stagnation line. 

These changes, first of all, only depend on the incident Mach num- 
ber. However, there is also a tendency that an increase in the 

angle of attack again leads to flow fields with a stagnation line, 
which is at least apparent for small Mach numbers. However, the 
proof that again leading edge vortices are created has not yet been 
given . 

For the flat pressure side, we can give the following physical 
explanation for the position changes of the stagnation line. For 
large Mach numbers, the shock is almost parallel to the underside 
in the central region and is only curved in the vicinity of the 

leading edge. In this outer region, the stream line has to be lo- 

cated whose continuation through the shock makes up the stagnation 
line. For a fixed angle of attack, the increase in the Mach num- /50 

ber then causes the shock to approach the bottom side. The almost 
plane central part of the shock is enlarged, whereas the curved 
part and therefore the stagnation line again are displaced in the 
direction of the leading edge. The reverse process can be observed 
both by reducing the Mach number or by increasing the angle of attack. 
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The position of the stagnation line can therefore be considered 
as an element which determines the lee side flow. The closer It Is 
to the leading edge, the less the underside chape can Influence the 
lee side flow. This Is clarified by a result from experiments. In 
the region of separation with Internal shock, a comparison of the 
flat lee sides showed that both the pressure distribution of the 
top side as well as the flow field (oil film photographs) are almost 
Identical for the delta model and the cone model. Here we can see 
that for higher Mach numbers, the stagnation line is Immediately 
adjacent to the leading edge and the influence of different under- 
sides (delta shape with V = 40° and conical with V = 75°) is not 
important for the lee side flow. 

Summarizing, we can say that depending on the incident flow 
conditions the shape of the bottom side has a non-negllgible In- 
fluence on the lee side flow. For the higher Mach number range, 
because of the position of the stagnation lines, the Influence is 
small, so that the top and bottom side can be considered separately. 


4.3.5. Additional remarks about the Influence of the cross-section 
shape 


We have previously discussed special geometric variables of the 
cross-section shape and in the following we will compare the Investi- 
gated shapes. 

An evaluation of the oil film photographs (according to Fig. 15) 
can be done here by interpreting the inclination angle of the wall 
stream lines in the vortex region as the intensity of circula- 
tion. With this assumption. Fig. 47 shows the vortex intensities 
on the top side of various delta wings as a function of the angle of 
attack. We find a good correlation of the values shown with the 
angle vp , which is measured in the cross-section plane between the 
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plane of the leading edges and the separation point. The angle 
Is used as an average value In Pig. 4? for the different 
lee sides and for a fixed angle of attack (« * 10°), even though 
in general It depends on the angle of attack because of the posi- 
tion change of the separation line. Nevertheless, up to average 
angles of attack (a ~ 12°) the vortex intensities of different lee 
sides can be compared. Higher angles of attack then lead to the 
already mentioned collapse of the top side flow and the values 
can no longer be correlated. 

In the previous sections, we have shown that the flow fields 
for example for delta shape lee sides basically do not change, but 
are only Influenced by means of the angle 0 , which changes the 
geometry. This is also supported by the’ separation line angle , 
which apparently is directly dependent on the shape angle for the 
plane top side shapes with 2-8 . Such a relationship can 

no longer be derived for curved lee sides, because then additional 
influences modify the flow field because of the shape. 

The effects of the different cross-section shapes of the delta 
wings can also be felt in a collective representation of the pres- 
sure distributions for 1 two angles of attack, shown in Pig. 48. In 
the separation region with leading edge vortex (M = 3.5 and a = 
5*5°), the pressure distributions are very different in the ex- 
ternal wing regions up to the leading edge, because the shape, 
position and intensity of the primary vortices vary greatly. In 
the central region, on the other hand, one finds attached flow for 
almost all wings and one finds comparable pressure coefficients. 

In the region of separation with internal shock (M^ =3-5, a = 

(i 

IQ ), the pressure distributions over the bodies with different 
angle 0 are very similar qualitatively, but do have clear quanti- 
tative differences. This has to be interpreted in the following 
way. In spite of different cross-section shapes, the individual 
elements of the flow field, such as expansion, internal shock, etc., 
are fixed more at certain positions in this region. The vortex is 
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no longer the dominant element and only has a small Influence on 
the pressure distribution. In contrast to this model, the delta 
wings have a different behavior with special cross-section shapes. 
The cone model with a stagnation point on the lee side only reaches 
the pressure level of the other wings In the central plane. For 
the modified cone model. In the vicinity of the leading edges we 
find pressure values which are comparable with other models. Be- 
hind the steep pressure Increase caused by the Internal shock, we 
find a similar tendency and a similar pressure drop as In the cone 
model. This indicates that for delta wings with cross-sections 
which are constructed from simple geometric shapes, in part the 
corresponding flow fields can be superimposed in order to obtain 
models of the flow processes. 

The aerodynamic coefficients are Important for the practical 
use of delta wings as aerodynamic bodies. As an example shown In 
Fig. 49, we will demonstrate how from the integration of the pres- 
sure distributions, certain lift coefficients of different wings 
will change with angle of attack. Four different cross-section 
shapes were selected and for two cases, we show the additional lift 
gain on the lee side. The comparison of the convex delta model and 
the flat delta model shows that in the region of separation with 
vortex formation (a < 5.5°), the lee side lift of the flat surface 
(0 ~ Q°) is greater than that of the delta-shaped wing (.0 = 9°). 

On the other hand, the values in this region of separation with 
internal shock are almost equal in size, as can already be con- 
cluded from the already mentioned similarity of the flow fields . 
Most of the lift coefficient is the result of the pressure side, 
of course, so that flat or slightly curved shapes will result in 
the highest value, whereas delta shape and conical undersides 
lie substantially below these values. Therefore, it is understand- 
able that out of the wings considered here, those with a slightly 
convex bottom side (0 = 5°) and a delta-shaped lee side (0 = 9°) 
will have the highest lift coefficient in the investigated angle of 
attack range. In order to evaluate the wing shapes, however, one 
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has to consider further aerodynamic coefficients and speed ranges. 
Therefore, it seems that this cross-section area is a favorable 
design for the underside and the top side with simultaneously an 
acceptable volume distribution. 

The lift coefficient which is also shown for two theories 
for flat delta wings '•shows good qualitative agreement for the 
linearized supersonic theory and good quantitative agreement for 
the "thin shock layer” theory. For a =0, we already have finite 
lift coefficients, because only in this way can one make a com- 
parison with thick delta wings where the underside already has an 
inclination towards the incident flow because of its form (0 > 0°). 

Additional information in Fig. 50 shows the different develop- /53 
ment of the vortex positions for different wing shapes . For the 
wave rider, the influence of the leading edge shock is very dominant 
for the lee side flow, so that compared with other wings, the vor- 
tex is displaced only slightly in the direction of the central plane 
for moderate angles of attack. Oil film photographs show that large 
area separations occur more easily than for other delta wings. On 
the other hand, there is a very stable flow field over the cone model. 

The internal shock changes its position hardly at all for average 
angles of attack, so that there is an almost constant vortex posi- 
tion. It is interesting that at the high angles of attack, in the 
present case a = 20°, the vortex separations are about the same for 
all of the wings. From this we have to conclude that here again 
we have a phenomenon which is independent of the delta wing geometry, 
which was already observed elsewhere. 

Summarizing, we can state the following important influences 
of cross-section shape on the lee side flow field: 

— For a fixed shape angle 0, the leading edge angle ^ has 
to be selected for stable flow states in such a manner 
that neither excessive flow angles occur at the front 
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side (total separation on the lee side) nor do negative 
values occur (stagnation point on the lee side). 

— Strong and discontinuous changes In the cross-section 
area can lead to substantial pressure Increases and 
temperature Increases. 

— By a favorable combination of the surface shape and the 
volume distribution, we find the maximum lift coefficient 
for the convex delta model (with delta-shaped lee side) 
out of the delta wings investigated here. 

The transfer of the results obtained here for practical applica- 
tions of thick delta wings as aircraft body forms has to be comple- 
mented by the following remark. Basically a certain volume has to 
be available for an aircraft. One characteristic quantity might 
be the volume parameter t sV/F 3/,! (V = volume, F = base area). 

If one assumes that in practice one requires at least a volume 
parameter of t = 0.08, and the values of t > 0.65 are no longer /54 

of Interest, then this can be specified for a wing with a delta- 
shaped cross-section by means of th«£ angle 0 . For a delta wing 
with a sweepback of K-l'f » we then find ~0'* 085 and 

T © = 20 ° =0 * 658 ' 

If we consider the example of a cruise Mach number of M = 4, 
then the flow region of interest on the lee side in the 
diagram can first be limited by curves corresponding to 0=8° 
and 0 = 20°, Fig. 51. In order to provide stable flow states with- 
out large area separations, we conclude from the present discussion 
that a maximum normal angle of attack of a Nmax = 30° cannot be ex- 
ceeded. In addition, experiments have shown that on the one hand 
no basic changes of the flow field can be expected when the angle 
© is varied and, on the other hand, for angles of attack of 
a ~ 25°, large separated regions occur on the lee side. In this 
case, there was a tendency to separations already for small angles 
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of attack In contrast to thin wings, so that there Is a further re- 
striction caused by the angle of attack ex v . As the lower limit 
for this example, we will assume the case where becomes nega- 
tive and a stagnation line with a high pressure is created on the 
lee side. Therefore, we have specified a Ninln * 0°. The remain- 
ing region inside the shaded lines is characterized by two differ- 
ent flow types, which are divided between separation with vortex 
formation and internal shock by means of the boundaries established 
in the present report. The sketched diagram is only an example for 
two parameters which have to be considered as restrictions for the 
top side in the design of delta wings in order to obtain stable 
flow states. 


5. Influencing of the lee side flow by disturbances in the /55 

tall region 

Some of the experiments showed that even for supersonic speeds 
the flow on the lee side can be influenced by disturbances down 
stream of the trailing edge. For example, experiments in the Cam- 
bridge wind tunnel at angles of attack of a = 15° showed an increase 
in the static pressure on the lee side, which sometimes was above 
the incident flow values [ 76 ] . Possible causes of this might be the 
interference or blocking in the tunnel or changes in the leo side 
flow, for example the bursting of the vortices and formation of 
shock systems. In order to clarify these questions and in order to 
avoid influences of the tunnel itself, we performed experiments in 
different wind tunnels. In the following we give a brief summary 
about these tests. 

As an introduction. Fig. 52 shows the pressure coefficient 
along the flat and delta-shaped lee side as a function of angle of 
attack. The results are from the three wind tunnels and show a 
pressure increase for certain angles of attack, which sometimes is 
above the values for the incident flow. One remarkable feature of 
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the results from the RAE tunnel Is that higher angles of attack 
again lead to underpressures on the lee side. A comparison of 
the suspensions of the models shows that In the Cambridge channel 
there Is a relatively thick sword with a wedge In the pressure 
side wake. On the other hand. In the RAE tunnel the wings were 
suspended on a wedge with a very large aperture angle, which 
penetrates Into the lee side wake as well as into the pressure 
side wake (Pig. 13)* The ILR measurements were carried out with 
a suspension configuration similar to that which was used In 
Cambridge and therefore show comparable results . Experiments 
for other longitudinal positions show that first in the vicinity 
of the trailing edge, the pressure values Increase but further 
upstream one still cannot find any influences. 

The static pressure on the underside of the wings, the pres- 
sure at the channel wall upstream and downstream of the model, as 
well as Schlieren photographs, prove that these phenomena cannot 
be contributed to a collapse of the flow in the test section. Ex- 
periments in different wind tunnels show that this cannot be ex- 
plained by specific tunnel disturbances. As Pig. 52 allows one 
to conclude, and from experiments with similar models 1411 , in 
which such a phenomenon was not observed, the different wing sus- 
pensions can be looked upon as a possible explanation for this. 
Pirst of all, we have to assume that there is an influence on the 

lee side flow because of disturbances from the tail region. 

In order to estimate possible disturbances of the lee side 
by the influences of the flow In the tail region, we carried out 
two groups of experiments: Pirst of all, the effects of perturb- 
ing bodies on the lee side flow were investigated and then we in- 

vestigated the influence on the tail flow. In the following, we 
give a brief summary about the results, and some of them were al- 
ready reported in [ 77 J • 

Pirst of all, experiments were performed with and without 
pressure hoses or coverings to simulate these hoses. We were able 
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to show that already by means of changes In the immediate vicinity 
of the tail, i.e., without any obstacles which penetrate beyond 
the top side, the influence on the lee side vortices is great. 

Similar observations were made by Richards [78] with a cone model 
(2). it was found that the vortices are non-conlcal for a certain 
angle of attack and then become conical again for higher angles 
of attack. Already Lukoslewicz [791 indicated conical separation 
on a Pitot tube. 

Further experiments were performed with wedges or blunt bodies 
[8Q] directly on the lee side or in the wake of the bottom side or 
the top side. Both the aperture angle and the height of the wedges 
as well as their position upstream and downstream of the trailing 
edge were varied. However, in these experiments it was not possible 
to obtain qualitative information about the influences of the dif- 
ferent parameters. In all of the cases with perturbing bodies, we 
found an increase in the static pressure as well as a bursting of 
the lee side vortex. The pressure Increase also occurred on one 
symmetry half of the wing even for asymmetric incident flow states, 
and overall it was dependent on whether the perturbing body was lo- 
cated in the wake of the underside or the top side. 

In addition, the effects of perturbing bodies on the tail flow 
were Investigated (see also [ 8l] ) . As a function of angle of attack 
in all cases there is a substantial increase in the tail pressure 
when wedges are installed either in the underside wake or the top 
side wake. 

6. Calculation of the flow field /57 


Up to the present there has been no theoretical calculation 
method known which gives a complete and realistic description of 
the flow field on the lee side in the region of "separation with 
shock." This is essentially based on the fact that very complex 
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processes occur on the top side, such as boundary layer separation 
with vortex formation and compression shocks with Interaction with 
the boundary layer, and this makes It difficult to obtain a closed 
system of equations which can be solved. 

It will be shown that a calculation based on a flow model de- 
termined from the experiments Is very difficult. For this reason 
It seems appropriate to Introduce models which are based on a 
highly simplified flow model. A comparison of the two calculations 
and a comparison with the experimental results will verify that the 
great effort to determine the lee side flow Is justified. 


6.1. Determination of the flow variables based on a highly- 
simplified lee side model 

In order to compare and estimate the present experimental 
and theoretical results, we will determine the Mach numbers, pres- 
sure coefficients and the position of the Internal shock using a 
simple calculation method. 

The flow model shown in Fig. 66 is similar to the one of 
Fowell [52], where In a plane perpendicular to the leading edge 
the flow variables are first determined. These are then the basis 
for the subsequent Prandtl-Meyer expansion around the leading edge. 
The additional component in the direction of the central plane 
which is created must then be equalized again by means of an in- 
ternal shock, as already mentioned. The Internal shock is assumed 
to be flat in this case and perpendicular to the wing surface. 

The Mach number perpendicular to the leading edge is already 
known as 

y 1 -sin 7 A.* cos 2 a 
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with the corresponding angle of attack 


a N * ore tan 


tan a 
cos A. 


With the assumption that the expansion Is complete, as already 
shown In Fig. 66, we obtain the Mach number M Ne from the Prandtl 
Meyer function 

# 
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The pressure coefficient Is then 
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The angle at which the flow Is Inclined In the direction towards 
the central plane Is found to be 

M Ne 


y s arc tan ( 
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Since this is also the deflection angle for the internal shock, we 
have 5 - v , and therefore from 


M e # y M J f ♦ ( M„ • cos a sin A )* 

we can determine the remaining quantities through the shock, such 
as the Mach number , the shock angle © and the static pressure 
ratio 
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The pressure coefficient is then 
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The position of the Internal shook is found from the angle which It 
forms with the symmetry plane Y-®-6 


and therefore 
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Prom these quantities, we have specified the pressure distrlbu 
tlon over the wing and this is then available for comparison with 
experiments and the calculation method now described. 


6.2. Determination of the flow field based on an experimental lee 
side model 


Before we discuss the developed calculation methods, we will 
first give the basic flow model. In order to remain within the 
region of separation with the shock and also to obtain comparison /59 
possibilities which are known from other calculation methods and 
experiments, the range of validity of the calculation Is restricted 
as follows. We will only consider delta wings with a plane flat 
lee side and sweep angles between A= 65 c and 85°. Over the entire 
flow region, we assume supersonic leading edges ( a =0°) so that we 
always have 

cos A. ^ , 

I 

Sm ^ 

where sin {*„ = l/M w 

The calculation method is therefore only valid for Mach numbers 
between = 3-5 and 7, as well as for angles of attack of a = 12° 
to 25°. The shape of the underside is then Introduced in the cal- 
culation only for determining the separated or the attached leading 
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edge shock, because In this Mach number range the lee side flow Is, 
for the most part. Independent of the underside (see Sec. 4.3.4) 
because of the fact that the stagnation lines are adjacent to the 
leading edge. The assumption of conical flow has been confirmed, 
for the most part, by experiments (see Sec. *1.1), so that in the 
calculation we can assume a cross-section area perpendicular to 
the leading edge. The experimentally determined flow model is 
shown in' Pig. 67 in this plane, and only for one-half of the wing 
because of symmetry, and it can be divided into the following 
regions . 

1. Incident flow in the cross section plane; 

2. separated leading edge shock; 

3. flow behind the shock perpendicular to the leading edge; 

4. boundary layer above the wing within which recirculation 
is formed; 

5. Prandtl-Meyer expansion over the boundary layer in the 
vicinity of the leading edge; 

6. internal shock; 

7. flow above the boundary layer which is thickened by the 
internal shock in the central region of the wing; 

8. parallel flow, at least in the symmetry plane. 

The flow field being described can be divided into two charac- 
teristic regions: The region of the boundary layer with friction 

is considered in the calculation as the surface of the wing, so 
that only the contour is determined, and the flow inside the bound- 
ary layer is ignored. When determining the Pitot pressure isobars 
in planes perpendicular to the wing surface, we introduce a simple 
trial solution for the state variables inside the boundary layer 
in order to have a possibility of making comparisons with the ex- 
periment in the region with friction. The frictionless region can /60 
again be divided into two regions, separated by a compression shock. 

The outer region of the wing is limited by the leading edge shock 
and the inner shock, whose shock front runs along conical lines 
through the tip of the model. A separation of the inner and outer 
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wing region exists because disturbances In the central region out- 
side of the boundary layer cannot propagate through the Inner shock. 
The calculation of the flow field Is done separately and In the 
order of regions shown above, Essentially, the theory contains 
the computation steps shown In the macro-flow diagram of Fig. 68, 
which will now be discussed. The Incident flow conditions and the 
geometry of the body are introduced as input data. After this, one 
decides about the leading edge type and the question of whether 
there Is an attached leading edge shock or a separated one. After 
the shock relationships have been used using the leading edge shock, 
the calculation of the expansion above the boundary layer has to be 
carried out in three dimensions. In order to finally obtain a veloc- 
ity field parallel to it in the symmetry plane. Since the expansion 
creates an additional velocity component, an internal shock is 
created from the mentioned symmetry conditions on the wing, which 
causes a deflection of the flow. The subsequent flow around the 
"thickened” boundary layer is calculated with different methods, 
depending on the speed range. In the supersonic range, the Prandtl- 
Meyer expansion is used, but in the subsonic range the panel method 
is used. The final flow field can only be determined by means of 
an iteration. 

Coordinate systems 

The coordinate system is selected according to various criteria. 
First of all, it should be possible to have a direct comparison be- 
tween theoretical and experimental results, but on the other hand 
one should obtain the simplest possible solution of the equations 
in the calculation method. 

First of all, we will select an aerodynamic and cartesian co- 
ordinate system whose origin is at the type of the model. The x- 
axis with a unit vector i Is parallel to the incident flow, the 
y- and z-axes with the unit vectors f and lc are shown in Fig, 69. 
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The; analytical treatment of the flow field is carried out in a /6l 

plane perpendicular to the leading edge of the model. Therefore, 
it seems appropriate to Introduce u body-fixed coordinate system, 
and the x-axls is parallel and the y-axis is perpendicular to the 
leading edge. This second coordinate system is obtained from the 
aerodynamic system in which the X-Z plane is rotated by the angle 
(90° -A.) and the X-Y plane is rotated by the angle a. The unit 
vectors here are X, m, n • The transformation equations from 
the body-fixed system to the aerodynamic system are therefore: 

X « x-cos a sin A.* y cos a cos A.-Z sin a 
Y * -x cos A.* y sin /V (1) 

Z « x sin a sin A.+ y sin a cos A.+ z cos a 

If we include the separated leading edge shock in the calcula- 
tion, we find a third coordinate system whose x-axls is parallel 
and its y-axls is perpendicular to the shock trace in the plane 
on the top side of the model. The rotation angles for a similar 
coordinate transformation as in (1) are the angle of attack a and 
the shock angle A-s • 




The leading edge shock 


The question of whether there is an attached or a separated 
leading edge shock at the leading edge is made based on the "thin 
shock layer" theory of Squire [82] . First of all, two new parame- 
ters are introduced by using a density ratio through a shock 
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i 

i 
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c = 


h 

e ,/2 • b 


O) 


From boundary layer considerations, we then find that in the 
C-fl plane, as shown in Fig. 70, the line Q, - C = 2 delimits sepa- 
rated and attached leading edge shocks. 


For the attached leading edge shock, we can have no influence 
of the top side on the bottom side, so that we can justify a sepa- 
rate calculation. The initial values for expansion at the leading /62 
edge are found directly from the values of the undisturbed flow. 

The Mach number components in the body- fixed coordinate system are 
obtained from scalar multiplication of 

M 0 , s m;n} , 

Therefore, the component parallel to the leading edge is 


M ot * cos a sin A. ( 5 ) 

and the component normal to the leading edge in the plane of the 
top side is 

M 0 m = cos a cos A ^ 

and the component normal to the top side is 


Moh^-Mc sin a 


(7) 


The Mach number component perpendicular to the leading edge 
is then the vector sum of ( 5 ) and ( 6 ) 


M hi = Mo,, • yl-sin 2 A cos 2 a 


( 8 ) 


The important angle of attack in this plane is also found from 
the Mach number components 
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tan a. 


Mqh w ton a 
Mqm cob ^ 


(9) 


For the case of a separated leading edge shock with - C) <2, 

we have to determine the shape and the position of the shock. The 

exact determination of the shock geometry is very complex even for 

simple bodies. For the present problem, we selected a simplified 

approximate solution. A first support of this are experiments with 

which the shock form was determined for different angles of attack. 

For a rhombic delta wing at M^ = 7, the results of Collis [48] are 

shown in Fig. 71. In the plane of the leading edges, the changes in 

the shock separation y/s = 1.1 at a = 0° and y/s * 1.2 at a * 17° 

or 25° are very small. Figure 72 shows an additional example of 

our own test results for M =3.5. In addition to evaluations of 

00 

Schlieren photographs, we also show results of the theory of Moeckel 
[833 and the "thin shock layer" theory. The latter has the disad- 
vantage that only values within the wing span y/s < 1 can be deter- 
mined and therefore the interesting region near the leading edge is 
not covered. On the other hand, the theory of Moeckel only gives 
the shock separation along the symmetry axis of the incident flow, /6 3 

whereas the shock shape is assumed to be approximately a hyperbola. 

Figure 73 gives an evaluation of this theory for two-dimensional 
and rotationally symmetric bodies in conjunction with our own experi- 
ments. It is exactly in the lower Mach number range that there is 
a strong change in the shock separation, whereas above M ~ 2.5 
there is approximately a constant value. If we use this knowledge 
and the experimental results, we can approximate the shock separa- 
tion referred to the span of the delta wing in the form of a hyper- 
bola by the following relationship: 

y S 0,748 . 

s ' M w -1 

Therefore, it becomes possible to calculate the shock trace angle 
A. s in the plane of the leading edges (see Fig. 69) as a function 
of the incident Mach number M^ and the angle of attack a. 
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(10) 


h.,.BK/2 » arc ton 


[' M.» 7 




cos 4 (la) -1 


1 ) tan I tc/2 -A.) 


The variation of this curve for different sweep angles Is 
shown in Fig. 7 1 * and is compared with available experiment data. 
Equation (10) shows that the influence of the angle of attack is 
small compared with that of Mach number which was already men- 
tioned in the discussion of Fig. 71. 


In addition to knowing the shock separation from the shock 
trace angle /\s , we only have to know the inclination angle of 

the shock with respect to the top side $ to continue the calcu- 
lation, which is shown in Fig. 75. This angle is determined in 
the plane perpendicular to the leading edge shock and is referred 
to the normal to the top side plane of the wing. The inclination 
angle in general differs only slightly from zero degrees, i.e., 
the shock is almost perpendicular to the plane of the top side. 

For the calculation, we set $ = 5° in all cases, which is an esti- 

5 

mat ion which seems to be justified based on experimental knowledge. 
This means we have specified all the parameters of the leading edge 
shock required for the calculation. 


In order to determine the flow variables above the separated 
leading edge shock, we define a coordinate system x, y, x in the 
plane of the top side in such, a manner that the x-axis is on the 
shock trace at an angle of A.*’ , as shown in Fig. 75. Just like 
in the attached shock, the Mach number component in the y-z-plane 
is found from 


Me; = M„ 


f 


sin 2 /u - cos 2 


(X 


(id 


and the effective angle of attack is 


_/6i 


cts 


arctan 


tan a t 
co* As 


(.12) 
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The shock angle for the use of the oblique relationships 
is 


Ps a 90*- Cf s ♦ $ s 


(13) 


Behind the shock in which it is to be considered plane, we 
obtain the following deflection angle 


6 _ sk /2 - arctan ton p ( 




M z 


£L 


2 (M|, sin 2 (3 s -t) 




( 14 ) 


and the Mach number 

ft 2 - 36 sin 2 ps-StMj sin 2 p s -1 )(7M| t sin 2 Ps + 5 ) ( 15 ) 

S ' «7M|, sin 2 p 5 ~l!(M 2 , sin 2 p $ +5) 

The total Mach number behind the shock in the X-Y-Z coordinate 
system is then 


M, = T(M w cos 2 a sir» 2 A,c,+ Nl s , cos (as + M 'os a cos A s 
♦ fvt s , sin (a s + Q s ) sin a ) 

-T( M„cos a sin A s cos A s -|Cl s , cos (a 5 +6 s ) sin A s ) 

♦ k( MojSin a cos d sin 2 /V s + £l s1 cos |a s +6 s ) sin a cos A 3 
- M s , sin ( a s -«- 5 S ) cos a ) 


( 16 ) 


In order to determine the initial quantities for the Prandtl- 
Meyer expansion, we make a scalar multiplication of Mach number 
M^ and the unit vector of the body fixed coordinate system. From 
the tangential component 

M|y = M, o t 

the component perpendicular to the leading edge in the top side 
plane 
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M lM = M,*>m 

and the component perpendicular to the top side 


we obtain from the scalar summation of and M^ N the magnitude 
M N2 the Mach number perpendicular to the leading edge. (If In /65 
the lower M*, range there is subsonic flow perpendicular to the 
VK, then nevertheless we can specify that the component required 
for the expansion Is at least = lj otherwise there would be a 
flow around the VK. Squire [10] also argued that because of the 
proximity of the VK shock, the sonic line starts at the VK.) The 
effective angle of attack is 

a N = urctan (Hm) , (.17) 

and the corresponding sweep angle is 

. / * Mjn \ ( "\ 

A.,. = arc cos ( — 

n ' CM,) ' 

In this way, from the oblique shock relationships we can determine 
all of the other flow variables, for example, Pitot pressure ratio 
and static pressure ratio. 

Boundary layer 

In the present calculation method, it will be found that it is 
not necessary to carry out detailed boundary layer calculations. 

As already mentioned, the displacement boundary layer is used as a 
fixed boundary for an isentropic (Prandtl-Meyer) expansion. Within 
the framework of this analysis, it seems sufficient to assume a 
simple two-dimensional boundary layer and to use the displacement 
thickness for a flat plate, as was done in l 84] 
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6 % £>( 0.664 * 1 . 73 ^) 

ifl U « /V aj 

Using the following rest variables, which prevail in the ILR tunnel, 
rest temperature T 0 = 293° K _ 

wall temperature If »T (1+0.169 Ml ) with t„ -t /(I + 4 — > 
speed of sound a,** 343.2 ~ (at T Q ) 

kinematic viscosity ,v OT 14,9 • icf 6 ~ (at T and P = 760 torr) 

i ° ° 

Sutherland constant r * l£_tL with S - 110° K) 

u T t-S 
w 

we obtain 


o* 

5 


:[0.GC4 




1,73 ( *,« -| 


■)/( 1» 0,1 G9 ML)] M i/2 V 


C • y » Vo 

2io 


where e = density ratio over a perpendicular shock. 


Prom experiments, we find that in addition to the Mach number, /66 
the angle of attack also has a substantial influence on the dis- 
placement thickness of the boundary layer. The ratio of these 
thicknesses for a wing with angle of attack and without angle of 
attack is found from 


6a__ 

65,0 


S 


The overall boundary layer displacement thickness is given by 

= 51 (2*43.3 a 9 ' 4 ) (19) 

In order to calculate the flow field, we require the contour 
of the boundary layer displacement thickness, which is first given 
by (19). 
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This relationship Is only used in order to determine the dis- 
placement thickness at a certain point. The position of the point 
can be found In the vicinity of the base point of the Internal 
shock and is described by means of the angle Y , which is between 
the symmetry plane and the corresponding trace line. This angle 
Is found from experiments based on an evaluation of oil film photo- 
graphs as follows : 


Y s * Uoo- A- ( 0.305 • a 4-0.083) (20) 

Here we have n® B ore sin (1/M^ ) , which is the Mach angle. From 

the angle y, we obtain the required distance from the leading edge 
Y so , which Is required to calculate the displacement thickness from 

y 50 ten ( 90-A -y) /p-i \ 

T" tan (BO- A 1 K 


This means that from equation (19), we can determine the dis- 
placement thickness ~ in the vicinity of the internal shock. 

Comparison calculations, however, give a negligibly small error if 
Instead of equation (19), one uses a simple parabolic shape of the 
boundary layer displacement thickness, using the coordinates at the 
point { y s * / 6* 0 ), determined above: 


6* s C s o f 


_y_i 1/5 

yso - 


(22) 


In this way we can determine the expansion from the displacement 
boundary layer determined here. 

Expansion in the leading edge region /&7 

The flow variables determined in the plane perpendicular to 
the leading edge are used as initial values for the expansion above 
the boundary layer. The parabolic displacement thickness along the 
span given by equation (22) is replaced by straight-line segments. 
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In steps, the flow variables in this region of the wing are deter- 
mined using Prandtl-Meyer expansion. 

The local flow angle <J> is equal to the angle for an attached 
VK shock for the first value of expansion. Also, the expansion Mach 
number here is M E1 * The corresponding position on the bound- 

ary layer Is determined using the inclination 

lfln, Pos B -57* i 6 c V<yy S o>' 5/6 (23) 

in such a manner that ip * \p os 

The division of the boundary layer into segments by specifying 
a fixed angular amount allows one to determine the velocity dis- 
tribution. The local flow direction is specified by 


<Pi « ‘PurA'P 

The Prandtl-Meyer angle is therefore 

V| b Vj. ( .+(*(>,- tp,.,) 


and the Prandtl-Meyer function is 


v = IT. arc tan arc tan 

and we can then determine the Mach number Mg^ 
as discussed in Appendix 1. 


V5*-i 

through an iteration 


This does specify the required flow variables in the cross- 
section plane, but we still have to introduce the tangential com- 
ponent into the calculation. This three-dimensional analysis is 
required because in the region of the symmetry plane, the flow di- 
rection must again be parallel to the Incident flow, i.e,, the scalar 
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multiplication of the velocity vector and the unit vector perpendicu- 
lar to the central plane has to equal zero. 

The vector of the Incident flow can be divided Into the follow- /6 8 
lng components 

U,* a Ut0 + UmO* UnO 
and also behind the leading edge shock 

If, “ Du * U m4 * Uni 

The tangential component Is constant over the shock and Is 

07, 

and we obtain 

u, a u w ♦ (u m i-iT ma i*(ir nl -u;, 0 ) (25) 

By dividing through the magnitude of the Incident flow speed, 
we find 

Gj_ _ Up, + ( Ufnl “Umol * (Uni ~Urq) 

lOLUuJ lull 11U 

(25a) 


The unit vectors are given by 


as well as 


--J22 s rn und ii m J- * m 

|U m0 | lUm.l 

und Mm s- *n 
lUnol I Uni | 


Por the y-direction, we refer these unit vectors m to the 
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absolute magnitude of 





lull 


' **«■» I **#■ I 

Umt # lUm lj 

|Umil |UwCi| 


JLJmiJ. 1 U ttiU ^ 

|Umol |U W | 


where 


|IT W |1 e Mni’CQS ^ . Oja. # Mst cos *P 
jOinol M sl ‘COS0f s o 5) Mst cos a s 


pf < 


(26) 


and 



cos cr 5 • cos A§ . 


Here, 4> is the local flow angle which is defined directly 
behind the shock with '•? ~ u 3 t 6., and which is reduced by A(J> 
during the Prandtl-Meyer expansion. Also, the Mach number M ^ 
directly behind the leading edge shock is given by equation (15), 
but during the expansion the prevailing Mach number M^ is sub- 
stituted in equation (26). 

If this trial solution is also Introduced In the z-dlrection /69 
and if we consider that 

JL. .7 

iu*i 

then for the vector in the outer region of the wing, we find 
-^2- = T* f m (g(M E )’ cos i? - cos Os) 

Iff- I (27) 

-fT ( g ( M r ) * sin c>-sin a s )j cos As 

According to equation (2 6 ) , the function g(Mg) Is defined by 
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In the aerodynamic coordinate system with the transformation 


and 


m * i cos a • cos A* ♦ j sin A$+ k sin a cos A.$ 
n s-T sin «*!T cos a 


we obtain the equation (27) shown completely In Appendix II. 


In order to determine the overall Mach number Mg in this 
region, we exploit the fact that the critical speed of sound does 
not change above the shock. Therefore, we have the following re- 
lationship between the local and the critical Mach number: 


where 


and 




12 


_ lu?! - !ik! . iHm! 
* V liu 



Mill* 


H-i 

2 



( 28 ) 


( 28 a) 


The internal compression shock 

During expansion over the boundary layer in the outer wing 
region, an additional velocity component is created in the direction 
of the central plane. Since symmetric flow fields occur for the 
delta wing under consideration, the direction of the velocity must 
be parallel to the incident flow at least in the symmetry plane. 


79 


Accordingly, in the central region of the wing, there has to be a 
deflection of the flow through a shock. The shock ^hape and posi- 
tion are not known a priori and are determined from an iteration. 

The different velocity ranges which are the result of the expan- 
sion (see Fig. 67) lead to a curved shock surface. Since we have /70 
assumed conical flow, the shock surface is formed by straight lines 
which all run through the model tip. In the following, we will de- 
termine the position and shape of the shock and the vector of the 
speed behind the shock. With the assumption that the tangential 
component over the shock is constant, we obtain the following for 
the speed: 


U 3 = U 2 ♦ ( 0 ,,-U.a) (29) 

where and are the components perpendicular to the shock 

front. In order to evaluate equation (.29), we first have to specify 
the direction of the normal components. In an arbitrary point on 
the shock surface, we can define a tangential plane using two vec- 
tors Ol and # which are perpendicular to one another. 

Oi = i cos y - } sin y 

& = i sin v sin 5 *j cosy sin $-k cos $ 

Here y is the projection of the shock trace angle in the x-y -plane 
and $ is the local inclination angle of the shock with respect to 
the y-z-plane rotated by y . The vector product is then the normal 
vector which is perpendicular to any point of the shock front 

d x & - % = i sin y cos $ + j cos y cos $ + 1< sin $ (30 ) 


The velocity components using equation (.29) are then found 


from 


and therefore 


U,3 = U, • t and Uc2 » U 2 • £ 


{U.j-U.aJSdJj-E- U 2 o*)o g 
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The Prandtl relationship gives the relationship between the 
normal components in front of and behind the shock 


*2 


If we again introduce unit vectors, then we obtain the follow- 
ing for equation (29) 


!L = J£ * (JU-iWrflVO’* 

jlU |IU , 2 2 'MUJ ' 


(31) 


With the relationship between the local and critical Mach num- /71 


ber 


M 






2 1 4 3izl j^ z 


and equation C27) 


we can then write 


|u 


Pi = a i + b j ♦ clT 


I u 


Ml . or. bT. ck. C^L_. (jto. - 8). e 


x If we use the condition that in the central region behind the 
shock the flow is parallel to the symmetry plane, i.e., 


Pi . « 
|UJ 


j = 0 


(32) 


we then obtain 




(.33) 


For the solution, we can express the term ( pL o^V in two 
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different equations, which when substituted In (33) then result 
In two equations for the unknowns M e2 , Y and $> . By using equa- 
tion (30, we have 

f~~j « ‘t? = a sin y cos $ + b cos y cos 4 ■* c sin $ ( 3*J ) 

I '•'(ID 


If we use the magnitudes of the Mach numbers, then we obtain 


ik 0 <• - ! Usd . . !l kL - Mfi . Uld (35) 

|U.| * Tu„l * I Ur. | M 2 -n 2 ‘ M 2 IUJ 


Then we obtain the following two equations 


2 M; a ■ 


b*M- 


xrl M e2 I U. 2 .I . cos Y cos «1> 
|U«| 


and 


(36) 


2 

- X4 V cos 2 y cos* $+ 7- cos y sin >1' cos i> 

2 M * 2 ~1 . b 

<• -- sin v cos y cos 2 $ 
b 


(37) 


By specifying an initial value, we can determine the unknowns 
in (36) and (37) using an iteration (Appendix III). The leading 
edge shock gives an additional velocity component in the negative 
y-direction, whereas the subsequent Prandtl-Meyer expansion gives 
components in the positive y-direction. Therefore, an expansion 
surface has to exist on which the direction vectors are parallel 
to the incident flow. Wherever the expansion surface intersects 
the shock surface, the shock intensity has to be zero. At this 
position, we have M e2 = 1 and u 2 0 f - 0, , so that from equations 

(3*0 and (35), we find 

a sin y cos $ + c sin $ = -L . /a 2 ♦ c 2 (.38) 

\A 2 ' 
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The velocity vector is 


- sp B °S0 > * C S0 ’ 


lu 


and Is parallel to the symmetry plane with the Inclination 


a so = arc tan (“$&) . 

' Q sc * 


The corresponding Mach number Mg forms a Mach cone, and the 
tip of the model or the coordinate origin is located on its genera- 
tor. In this way the shock angle y SQ can be given as follows in 
the x-y-plane 


Y SO c arctQii y tan^iiso -tan^ 

The inclination angle $ has to be determined from a solution 

5 O 

of equation (38) . If we calculate the corresponding positions in 
the coordinate system, then all initial values are given for itera- 
tion to determine the internal shock. 

In this iteration, we first assume that the inclination of the 
shock $ is constant, or the value is used from the previous cal- 
culation. From equations (36) and (37), we then find an addi- 
tional point of the shock surface in the cross-section plane. The 
connection line between the new and the old shock points produces 
the inclination angle which has to be varied until the value used 
in the calculation equations and which is the result of the geometry 
coincides within certain limits. 

The shock shape and positions can then be determined in the /73 

manner described. The velocity vector behind the internal shock 
is then found from equations (.31), (35) and (.36) 
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JtLi s qT+ b7* clT •» r ( i sin v cos $ ♦ j cos y cos $ + iT sin $ ) 

|U J cos \ cos $ 


from which it follows that 


j~-j = •' (Q-b tanyJ+k (c-b 


tan $ . 
cos y 


( 39 ) 


The magnitude of Mach number is determined in a similar way 
as was done with equation (28), so we obtain 


where 


M, 



i 1 ' 2 - m: (i 


*K-1 

X*1 



(40) 


M,* * 


IU3I . luj . ijjcoi 
a* IUJ a* 


and |U C0 j /a* is defined by equation (, 28 a). 


The pressure increase above the internal shock is the reason 
for the thickening of the boundary layer in the central region. 
Here there is an expansion which results in a velocity component 
in the direction towards the symmetry plane. Therefore, the equa- 
tions for the internal shock have to be expanded so that a larger 
deflection of the flow is achieved and so that only occurs 

after the expansion. 


Using the condition U- D T^ 0 ""^ » Instead of equation (33) 

we obtain 


and therefore for equation (. 36 ), we obtain 


( 41 ) 


2 _ tb«K) M; 

■H4l M e2 IU 2 I . cos Y cos $ 

1 '-'Oil 
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and finally the velocity vector behind the Internal shock is 


= I (a - (b+ K! tan v j i b -( b+K J) 


♦ k ( c - ( b + K ) 


tan 4~ 

cos y 


( 42 ) /n 


With an appropriate Iteration, we have to determine the value 
K in such a manner that after the internal expansion, one obtains 
a vector parallel to the symmetry plane, as discussed above. 


The inner expansion 

The displacement boundary layer downstream of the internal 
shock is again considered as a fixed body surface at which the in- 
ternal expansion occurs. The form of the boundary layer in this 
region is given in various ways (.Cross, Collis : Pig. 7) and will 

be described by means of an ellipse with a principal axis in the 
symmetry plane. The equation of the ellipse can be determined from 
two points. Prom the intersection point between the internal shock 
and the boundary layer, we find the "separation point" which is 
also on the thickened boundary layer. The second point is defined 
by the displacement thickness in the symmetry plane &cl which 
was determined in various experiments. Figure 76 shows the de- 
pendence of the thickening on angle of attack for different Mach 
numbers. As an approximation, we can obtain the following equation 
from this : 

6 * 

« 0,97 • M^-a ♦ 0.001 M* ( 43 ) 

If the coordinates of the separation point S( y 05 /z 0S ) 
and of the point on the symmetry plane S CL ( y CL / 5^} are known, 
then the ellipse equation is 
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where 


C 2 


4 


(43a) 



1 


o* = t y os - y ct . ) 2 * 


Ca 


6a" 2 2 


os 


) 


In the continued calculation in the central region of the 
wing, we have to distinguish whether the velocity components 
reach supersonic or subsonic speeds in the cross-section plane. 

We will discuss both cases. 

In the case of supersonic speeds, the calculation is similar 
as already given for the expansion above the boundary layer in 
the outer region, in which first the ellipse is replaced by small 
straight-line segments. Starting from variable^ behind the in- 
ternal shock, the Prandtl-Meyer expansion then gives the Mach num- 
ber Mjj e and the local flow angle . 

The velocity vector and the total Mach number are then found 
from the following trial solution. The velocity vector behind the 
internal shock was 

iL.T‘d*k.e 

|U«,| 

where d and e are defined by equation (39)* 

In the cross-section plane perpendicular to the inner shock, 
we than find 


ri t (d* cos a sin y •* s sinasiny) 

♦ m t {d cos a ccs y * c sin a cos y ) 

♦ 7T t {-d sin a < a cos a ! 
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( 44 ) 


ORIGINAL PAGE U 
Of POOR QUALITY 


Since again we assume that the tangential component above 
the Internal shock does not change, we obtain the following for 
the speed 


U**U,*(U 4 t -U 1# ) 


If we define 



m 4 cos tp 4 -n t sir, tp 4 


and If U 3e/ |U M f is given by the m^- and n^ -components in equa- 
tion (44), then we find 


with 


IUJ 


f- I m 4 cos <p 4 -ivsin «p 4 )“S(M 4e ) 

- rn (d cos a cos y + e sin a cos y) 
-"n !c r.os a *• d sin a) 


n ( m 4# j 



In this way, we can write the velocity vector in the central 
region of the wing as follows : 

Jli. = f .{ ri cos a sin v + c 5,n a siM ^ * 

|IU ‘ 

,tn t !g! M 4i ) cos ip 4 J-n lg(M tc )*sin q> 4 ) (45) 

The complete equation in the aerodynamic coordinate system 
is contained in Appendix IV. 


The magnitude of the Mach number is found from 


1 


M, 


(*) 


1/2 


m; n 


% M r> 


#2 i - 1/2 


(46) 


with 


M* 


|UJ - lu 


Q 


1 ^4 I 

IUJ 


IU Q 


For low Incident speeds, around ~ 4, we obtain subsonic 
speed In the central part of the wing In the plane perpendicular 
and downstream of the Internal shock. The calculation of the flow 
variables in this case is done using the panel method and will now 
be described. 


The velocity field can be represented as the sum of the veloc- 

-V -V 

lty behind the shock U, and the perturbation velocity U given by 
the elliptical thickening of the displacement boundary layer. 




u 3 4 U s 


The calculation of the perturbation field U is done using 
potential theory, and a disturbance potential is introduced, which 
for example can be produced by a source distribution on the surface 
of the boundary layer. Then we can establish an integral equation 




W) 


whose physical meaning is shown in Fig. 77. The right side of this 
equation is the normal component of the incident flow at the point i, 
The left side, first of all, consists of the term 2ru\. , which 

is the normal velocity component directed outwards, caused by the 
source intensity in the immediate vicinity of the same point. The 
surface integral, on the other hand, is the component which results 
from the source distribution at the points j on the remaining body 
surface at the point i. 
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Since equation (47) cannot be solved analytically for arbi- 
trary body shapes, the panel method [851 Is used, and the body 
surface Is divided Into a number of plane surfaces. On each of 
these panels, one assumes a constant source Intensity with an un- 
known value, and then q^ can be placed ahead of the integral. The 
integral over the Individual point Is then reduced to a geometric 
problem which can be solved simply. We then flxid a system of 
linear equations 


N 

^ v\i 

i=l 


n: • U, 


(48) 


ni 


which is an approximation to the Integral equation (47). The 
physical meaning of the Individual terms is also shown in Pig. 77. 
An additional description of the method used here is given in [861. 


The solution of equation (48) gives the velocity field in the 

cross-section plane so that together with the tangential speed, the 
*+* 

total vector can be represented by 


j™, = *4 (c5 cos a sin y + c sin <x sin y) 
^ m A ( U s 0 ir » 4 ) ^ n 4 ( U s ° rf ; ) 


(49) 


The corresponding Mach number distribution is then found again from 
equation (46). 


In the previous section, we already discussed that the flow 

around the boundary layer thickening induces an additional component 

«»♦» **■ 

in the direction of the symmetry plane and therefore U, 0 j / 0 ist. 

The shock has to be iterated in such a manner that it produces an 
additional component in the direction of the leading edge. Accord- 
ing to equation (42), in the calculation of the inner expansion 
we have to include the j -component of the vector U^. The speed 
Ujj then only changes in the direction. This equation is also 
discussed for the aerodynamic coordinate system in Appendix V. 
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Static pressure and total pressure in the flow regions 


In the previous sections, we discussed the velocity vectors 
and Mach numbers In the entire flow field, and therefore we can 
now calculate the static pressures and the total pressures In the 
Individual regions. In the following the Index («>, 1, 2,3,4) re- 
fers to the flow region and (p) refers to the Pitot pressure behind 
a normal shock. 

Behind a leading edge shock , we obtain the following for the /78 
static pressure ratio: 

P, e 2x M s rsin* (50 \ 


For the total pressures, it makes more sense to give the ratio 
of the Pitot pressures in order to make comparison with experiments 

p p» .. ZlI , H . Jk (51) 

5“ ' p . p ~ p >- 

so that according to the shock relationships, we have 

?r’i 

n 

P L * f (Mji ' sin 0 S ) 

Pa* 

Pp 

— ■MM*,) 

In the outer expansion region in the vicinity of the leading 
edge, the static pressure ratio is found from a series expansion 
for an expansion over small deflection angles Av . The Mach num- 
ber is the one of the preceding expansion 
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from which we find 
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For the Pitot ratio, we can then essentially use equation 
(51). The outer expansion only causes the Mach number behind the 
leading edge shock to change and therefore we have 
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where 


f?'L 4 f ( M p i 
Fa 

and P 2 / p x = 1 Is an isentropic state change. 

The static pressure ratio over the Internal shock Is again /79 
given by the shock relationships and Is the following, referred to 
the incident flow pressure: 


El 2jl . 2l 

Poe *P 2 Pro 


(5*0 


where f ( /v W 3 where M e g is the perpendicular component 

ahead of the internal shock, The Pitot pressure ratio can then be 
determined from 


Ppi. „ fr . h, . JL . fa 

Pp ro " TFj P 2 Ppi PPeo 


(55) 


where 
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and 


!r- * hm„) 

ri 

^ * 1 (M t: ) 


where M e3 Is the normal component behind the Internal shock. 

In the region of Internal expansion in the vicinity of the 
symmetry plane, as already discussed In the previous section, 
there will be subsonic speeds and supersonic speeds. Therefore, 
In the case M^ e > 1, we use the equations of Prandtl-Meyer 
expansion for the static pressure, as already done in (52). 


If subsonic speed prevails Mjj <1, then we find 


n 






(56) 


The Pitot pressure ratio In supersonic flow is 


pa „ ^pa 
P p ' P 4 

* 09 


P, P„ 
Pa Ppa Pp ( 


^A ’ i ‘Pi 


where 


!M, t ) 

Pa 

and Pjj/p^ = ;1, because here again an lsentropic state change has 
occurred. For the case where subsonic Mach numbers Mjj e < 1 exists, 
the Pitot press. re ratio does not change and then corresponds to 

1 'p3 //p P w ' 


In all regions of the wing, we wrote the static pressure as a 

coefficient C to better compare it with experiments, which is 
Jr 


/ 80 * 
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defined by 


Cp 



0.7 


This means we have now given the most important flow vari- 
ables for a comparison with experimental results and other 
theoretical results on the lee side. 


7 . Discussion of the theoretical results and comparison with 
experiments 

In the discussion of the theoretical results, we will first 
describe the influence of the simplifications used in the calcula- 
tions in the form of empirical equations or fixed data. For an 
example, we will then describe the complete theoretical flow field 
which was determined and will carry out a variation of the most im- 
portant input variables. Finally, we will give a comparison with 
other theories and experiments, including those which are outside 
of the definition region of the present calculation. 

The empirically determined equations which were used in the 
calculation method agree quite well with the various experiments, 
but at this point we have to examine whether deviations have an 
effect on the theoretically calculated flow field. Four components 
of the flow model have been determined empirically: thickening of 

the boundary layer behind the internal shock, displacement thick- 
ness ahead of the internal shock, inclination angle <i> and trace 

s 

angle of the leading edge shock. 

The thickening of the boundary layer downstream of the internal 
shock in the vicinity of the symmetry plane has a very small effect 
on the remaining flow field. We will discuss the effects on the 
pressure coefficient due to this when we make comparisons with ex- 
periments . 
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The determination cf the displacement boundary layer upstream 
of the internal shock was based only in part on empirical knowl- 
edge, because the thickness was derived from the simplified bound- 
ary layer equations [84]. This trial solution seems to be justi- 
fied for the flow model, as the evaluation of control calculations 
in Pig. 78 shows. Here, we summarize the influences of various 
displacement thicknesses 5 S * 0 / s on the pressure coefficients and 

the position of the Internal shock. In addition to the value 6* 0 /s 
from equation (19), which in the figure is indicated by (j) , we 
show a comparison value which is about twice as large. Indicated 
by © . Prom the sketch on the left of the figure in the upper 

diagram, we can see that the initial point of the internal shock A 
varies greatly in the z-dlrection, which however barely has an in- 
fluence on the further calculation. The base point of the internal 
shock F (a measure for the trace angle y) changes only slightly in 
the y-dlrection. In the upper diagram, we show a pressure coeffi- 
cient upstream of the internal shock which changes only a maximum 
of 3 % with respect to the value at (j) . The greatest influence 
when varying the displacement thickness is exerted on the pressure 
coefficient just downstream of the internal shock. Here we can 
observe deviations from the value at (f) which amount to between 
5 % and 1Q$, which however still lie within the order of magnitude 
of the measurement accuracy of the experiments. 

The influence of the inclination angle $ of the leading edge 
shock is treated in the same way in Pig. 79. The average empirical 
value is $ s . « 5°, which was used in the calculations, and it is 
marked in the figure. The lower diagram shows the effects of vari- 
ous inclination angles on the position of the internal shock. The 
initial point A changes its position in the z-direction only incon- 
sequentially. Also, the base point of the internal shock is almost 
constant and only at large inclination angles $_ can we find devia- 
tions of around 10'# in the average value. This is also true for 
pressure coefficients in front of and behind the internal shock 
which are shown in the upper diagram. As already shown in the dis- 
placement thickness variation discussion, influences on the pressure 


coefficient are greater downstream than upstream of the internal 
shock. Nevertheless, we can still find a small deviation in the 
pressure coefficients compared with the comparison values. The 
angle 4> g ■ 5° can therefore be further used in the calculation 
equations, because even for changes in the real flow by A<S> s » + 5° 
we only find small influences on the flow field determined from 
theory . 

The greatest effects on the pressure coefficients and posi- 
ticm of the Internal shock are obtained when one varies the trace 
angle A, s of the leading edge shock. This is understandable if 
we consider that the values perpendicular to the shock plane form /82 
the departure values for the entire calculation. Figure 80 shows 
an evaluation made for the other parameters, and the trace angle 
A. s was varied between 65° and 70°. The value determined from 
equation (10) for the prevailing incident flow conditions, and 
geometric conditions is /V s * 67.3° . From Fig. 73, in which accord- 

ing to equation (10) the angle a. s is represented as a function 
of Mach number, we can see that for the selected delta wing sweep- 
back, in the range a. s > 63° .• small changes in the trace angle lead 
to a strong increase in Mach number. For i\ s >63° , a change in 

the trace angle has a small influence. This result is also re- 
flected in the lower diagram of Fig. 80 for the position of the in- 
ternal shock. For A. s >£8° , the initial point and base point change 
more with increasing trace angle than in the range A. s > 68° . For 
the pressure coefficients Cupper diagram), this means that the values 
ahead of the internal shock are influenced relatively little, but 
according to the position change, strong pressure differences down- 
stream of the internal shock occur. This means that we have a sub- 
stantial influencing of the flow field by the parameter A. s , so 
that it seems unavoidable to carry out the most accurate calculation 
of this value possible. On the other hand. Fig. 73 has shown that 
with equation (.10), we can obtain quite a good approximation, which 
results in a deviation of AA. s s= + i’ with respect to experiments. 

The resulting error is Justifiable for the pressure distributions 




but for the base point of the Internal shock, we obtain position 
changes of a maximum of 30# referred to the comparison value. 

At this point, with the available Information we cannot yet decide 
about the feasibility of the calculation method. Only a subse- 
quent comparison of theoretical and experimental results will allow 
a complete evaluation. 

The theoretically determined flow field on the lee side of a 
delta wing Is shown in Pig. 81 for M w = 3-5 and a = 14. 5°. For 
the selected incident flow conditions, we have a separated lead- 
ing edge shock, through which the velocity component In the cross- 
section plane undergoes an additional deflection. Along the dis- 
placement boundary layer, the flow expands, and over the wing, 
regions with different velocity result. On the first expansion 
wave shown, we find the initial point A of the Internal shock, 
which here at the same time is part of the Mach cone of the inci- 
dent flow. The further course of the internal shock is described 
by a curved line. Downstream of the internal shock, we again have 
a deflection of the flow because of the thickened boundary layer, 
until symmetry conditions are achieved in the central plane. /83 
Arrows indicate the most interesting positions and they indicate 
the size and direction of the speed. This flow field is, in princi- 
ple, maintained for all of the incident flow conditions specified. 

Figure 82 shows the pressure coefficients over the span from 
experiment and theory at two Mach numbers. We can see that for an 
average angle of attack (a „ 12°), the pressure coefficients are 
comparable both in the expansion region and in the central region. 

At high angles of attack (a ~ 20°) in the expansion of the wing, 
we still find a similar pressure level for the experiment and the 
theory, but from the theory we find pressure increases over the in- 
ternal shock which lire too high. In contrast to the theoretical 
results, the pressure distributions which are the result of experi- 
ments are continuous over the span. Discontinuous pressure in- 
creases from the calculation are obtained in the region of the 
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internal shock. In the real flow, there is a pressure equalization 
inside the boundary layer In this region. If we also assume that 
the shock in the vicinity of the boundary layer consists of indivi- 
dual compression waves, then we find the variation of an isentroplc 
compression which is Indicated in the upper diagram. 

When introducing the calculation method in the previous chapter, 
we first presented a simple method based on expansion equations and 
shock equations. The flow model was greatly simplified. For M ro = 4 
and a = 12.2° we then obtained a result which is also given as an 
example in Fig. 82. It will be shown that compared with experiment, 
we have an expansion which is too great and that the shock positions 
and pressure Increases are not comparable. 

A comparison with the theory of Babayev [53] is shown in Fig. 

83 . On the basis of the inclusion of the displacement boundary 
layer in the previous report and the centered Prandtl-Meyer expan- 
sion at the leading edge in the theory of Babayev, we find substan- 
tial differences in this wing region. Ahead of the internal shock, 
we achieve almost the same pressure coefficients but different posi- 
tions and intensities of the internal shock. The other pressure 
distributions shown are quite close to one another in the expansion 
region and are very different in the central region. On the other /8*l 
hand, the linear theory is already outside of its validity range 
(M^ < 10) and results in pressure coefficients which are below those 
of a vacuum. Oil film photographs of Rao and Whitehead [5] (with 
somewhat changed sweepback and angle of attack) give an indication 
of the position of the internal shock. This is indicated by an 
arrow in the diagram. We have shown, therefore, that the individual 
theories- give completely similar results. However, in this calcula- 
tion example, we have no experimentally determined presssure dis- 
tributions and we cannot evaluate the individual trial solutions. 

In the previously discussed pressure distributions, it has be- 
come clear that the vortex position and shock position cannot be 
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determined uniquely. The pressure distributions show a pressure 
equalization over the boundary layer and the oil film photographs 
also only give wall stream lines within the boundary layer. It 
can be assumed that the position of the vortex and the shock will 
occur In the same direction when the Incident flow conditions are 
changed. On this basis. Pig. 84 gives a comparison of the posi- 
tion of the base point of the Internal shock and the positions of 
the vortex centers determined from oil film photographs. The 
basis for this are the results from the CUED wind tunnel shown in 
Pig. 32 for a fixed Mach number and different angles of attack. 

We find a good agreement between the change In the vortex center 
and the internal shock in the span direction. The assumption of 
Cross [^ 9 ] that the separation lines are displaced In the 

direction of the leading edge Is therefore supported. This is not 
necessarily contradictory to the available pressure distributions 
If one assumes that because of the thicker boundary layers, the in- 
fluence of a shock on the surface is greatly decreased. The pres- 
sure Increase in the central region Is created essentially by the 
vortices inside the boundary layer and is substantially below the 
values which were achieved above the internal shock. We have to 
point out that with increasing angle of attack, the shock Intensity 
In the present calculation as well as according to the theory and 
experiments of Powell [52;] increase. 

When one varies the angle of attack, one finds quite good 
agreement (see lower diagram In Pig. 84) between the height posi- 
tions of the vortex center and of the base point of the internal 
shock. Based on the flow model, it is understandable that the 
base point of the shock will always be above the vortex center. 

Prom the previous results, we can already see that the calcu- /85 
lation results in usable results even for angles of attack of 
a < 12°, i.e., outside of the definition range. As an example of 
this, we show the experimental Pitot isobars determined by Monnerie 
and Werle [75] as well as the theoretical flow field. Fig. 85 . As 
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already described above with the ILR results, even here the Pitot 
pressure values show a separation bubble (Manrterie and Werle, 
separated region) at the leading edge and attached flow in the 
central region, and their positions from oil film photographs and 
Schlleren photographs are also indicated. Prom the theory, we 
find that the displacement thickness ahead of the internal shock 
is almost identical with the dimensions of the separation bubble, 
but downstream of the Internal shock we assumed a boundary layer 
(i.e., a thickening) which is parallel to the surface. Overall, 
we find good agreement between the expansion region and the posi- 
tion of the Internal shock, but for the central region we cannot 
make any determination except the indication flow direction (In a 
plane perpendicular to the internal shock) . 

Finally, in Fig. 86 we show how the flow field changes when 
one varies the angle of attack or the incident Mach number. For 
a fixed angle of attack = 14.5°, the figure shows on the left 

side the Mach number which was increased from Moo = ^ to 7. The 
leading edge shock moves in the direction towards the wing and 
therefore slightly increases the normal angle of attack (Aa^ ~ 3°). 
The simultaneous increase in the normal component of Mach number 
brings about a strong increase in the deflection angle through 
the leading edge shock, so that a stronger expansion is necessary 
until that expansion surfs, ce is reached on which the internal 
shock has zero intensity. The expansion of the internal shock in 
the z«dlrection is smaller based on the thicker displacement bound- 
ary layer, whereas at the same time there is a slight displacement 
in the direction of the symmetry plane. 

As already mentioned, the influence of the angle of attack on 
the leading edge shock is small, so that for it we have an almost 
constant position on the right sketch of Fig. 86. The variation 
in the angle of attack from a = 12.5° to 20° brings about a rela- 
tively continuous increase in the normal angle of attack, whereas 
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the corresponding Mach number changes only slightly. Based on 
these conditions, even behind the leading edge shock one obtains 
an almost continuously changing flow parameter collection, so 
that the Internal shock is displaced uniformly in the direction 
of the leading edge. In this way we were able to show how the 
flow field behaves when one varies, two Incident flow parameters. 

The trends which were found have been confirmed by experiments. 

Summarizing, we can say here that the empirically determined 
equations used in the calculation for the most part have a small 
Influence on the most Important quantities of the flow field. One 
exception to this is the trace angle of the leading edge shock , 
which greatly influences the position of the internal shock. How- 
ever, comparisons with the experiments have shown that the position 
of the internal shock is well represented by the present calculation 
method. As far as the pressure distribution is concerned, over the 
entire definition region we find good agreement for the pressure 
level in the expansion region. The pressure values downstream of 
the internal shock in the central region of the wing are comparable 
up to a maximum angle of attack of a ~ 15 0 . The relatively high 
complexity of the calculation method compared with other theories 
with greatly simplified flow models is therefore justified, because 
both the flow field as well as the pressure distribution, for the 
most part, agree better with experiments. 


8 . Summary 

In the present report we discussed lee side f* w over delta 
wings under supersonic conditions, both experimentally and theoretic 
ally. Since not much work has been done in this area, it was the 
purpose of this investigation to review present-day knowledge and 
to extend the information known about the flow field. 


The experimental part included thick delta wings with different 


r 




cross-section shapes, and the experiments were essentially carried 
out in two different wind tunnels . The experiments with the delta 
wing with a flat top side were used to analyze the flow types for 
the different incident flow conditions . The Stanbrook-Squlre 
region was used as a basis of this, which separates leading edge 
vortex formation and shock-induced vortex formation. The evalua- 
tion of the experiments showed the following: 

0 Confirmation of the region to the left of the Stanbrook- /87 

Squire region (M m = const, a * variable) with leading edge 
separation into primary and secondary vortices. However, 
at the higher angles of attack, we found a change in the 
type of flow. Here we find a shock which is almost 
parallel to the surface in the symmetry plane, which is 
the result of the downwards directed velocity components 
of the primary vortex. Instead of the secondary vortex, 
we only observe a dead-water region here. 

0 Inside the Stanbrook-Squire region CM„ = variable, a = const) 
the flow field changes continuously, so that no fixed bound- 
ary for the change in the flow types can be given. However, 
the most noticeable change in the experiments was observed 
around the center of the transition region (JV^ » 3-0). 

0 To the right of the Stanbrook-Squire region (H» = const, « = 
variable), we again found leading edge vortex formation for 
small angles of attack. When the angle of attack is in- 
creased, we find a flow type which deviates from the shock- 
induced vortex formation flow type after crossing a transi- 
tion region. Here we find a separation bubble at the lead- 
ing edge and an Internal shock on top, for the most part in- 
dependent of one another. This type of flow is called 
"separation with shock" and lies to the right of this 
Stanbrook-Squire boundary only for average angles of attack, 
because at higher angles of attack we again assume vortex 
formation from the leading edge. 
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The experimental Investigations with delta wings for non- 
flat (conical, delta-shaped) top sides were performed In the 
"separation bubble with Internal shock" region. We were able 
to establish the following: 

0 The lee side shape gives changed Incident flow conditions 
In the cross-section plane. From this we find attached 
flow at the leading edge for the delta-shaped lee side 
together with a supersonic expansion. Downstream, the 
internal shock induces a separation with vortex forma- 
tion. 

0 The angle 0 between the ridge line and the plane of the 
leading edges can be defined for flat top sides and can 
define how "greatly delta-shaped" a lee side is. A change 
In this parameter only has a small Influence on the geome- 
try of the flow field. 

0 In general the cross-section shape of the delta wings in- /88 

fluences the aerodynamic properties very greatly. Excessive- 
ly large angles ¥ between the top side and the bottom side 
perpendicular to the leading edge lead to stagnation points 
on the lee side, whereas the flow leads to complete separa- 
tion over the wing for small angles V . 

0 For separated leading edge shocks, the shape of the under- 
side has an influence on the lee side flow which should not 
be ignored. However, in the higher Mach number range, be- 
cause of the position of the stagnation lines on the under- 
side, it is justified to treat both sides separately. 

In a few experiments, we observed strong pressure increases 
above the entire span which had not been previously observed and 
therefore we made additional investigations of this phenomenon. 
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For certain Incident flow conditions, we found that even for super- 
sonic speeds, disturbances In the tall region can influence the lee 
side flow downstream of the trailing edge. 

The theoretical results were restricted to the "separation 
bubble with Internal shock" region and are based on a flow model 
which was determined from the available experimental data. Em- 
pirically determined equations for the position and shape of the 
leading edge shock, as well as for the displacement boundary layer, 
were substituted. The resulting calculation equations for the 
velocity distribution over the lee side had to be solved using an 
Iteration process. Comparisons with experimental pressure dis- 
tributions and flow fields show good agreement for the expansion 
region and for the position of the Internal shock up to medium- 
sized angles of attack (a ~ 15°)* At higher angles of attack, on 
the other hand, the pressure jump over the Internal shock is very 
strong and can hardly be compared with the measured pressure dis- 
tributions over the span which are almost uniform. However, the 
flow conditions are, for the most part, not clarified for this 
region. 


For future work, therefore, we have shown that the Stanbrook- 
Squlre region Is not the only limit between the various flow types . 
In particular, we should point out the region of small or very large 
angles of attack, in which apparently substantial changes in the 
flow field exist. In addition, we have to point out the influence 
of the Re number, which in the present report was only discussed 
with regard to the different wind tunnels. In addition, the knowl- 
edge of the flow fields allows one to investigate the Influence of 
the tail region disturbances on the lee side flow. We have not 
obtained a thorough understanding aboxit the mechanism of this 
phenomenon. 
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APPENDIX 1 /101 


The Iteration of Mach number from the Prandtl-Meyer function 
for specified angle v first requires the exact knowledge of the 
function v = f (M) . The first and second derivative is given in 
the following for * = 1.4 : 

v s "fs- arc tan -arc tan / M 2 -T 



dM 2 M 2 1 M 2 + 5 } 2 -(M 2 -1) 


j 



From this it results that for v = 0 the function runs tan- 
gentially into the M-axis at M = 1 and has an Inflection point 
for M w = 1.5812 (at v ) . Therefore, using the following flow 
diagram, any Mach number can be determined from the specified 
Prandtl-Meyer angle 









— . 1 

Here £ Is a freely selectable accuracy constant. 

APPENDIX II /IQ 2 

The determination equation for the vector In the outer ex- 
pansion range for the aerodynamic coordinate system according 
to equation (27) Is: 

U 2 T* r 

jy“j = ' [Ikos a cos A cos A. s (glM e )‘Cos ip-cos ct s ) 

♦ sin a cos A. s (g (M E )'sin ip -sin ct s )] 

♦j [sin A cos A. 5 (g(M E ) cos ip - cos a s )] 

♦k [sin a cos A. cos A. s (g(M E ) cos ip - cos a s ) 

-cos a cos A. & (g (M E )-sin ip -sin « s ) j 

APPENDIX III 


Summarizing equations (36) and (37) results in an expression 
which only depends on the shock parameters Y and 

cos 2 $(cos 3 y+ — r sin 2 y cos y)+cos 2 y(2 ~r~sin ycos 2 $ 
b‘ b 


+ 2 ~ cos $ sin $) 1 2 ™ sin y cos y sin $ cos $ - — sin y 

- cos y ( ~~ + — sin 2 $♦ IML U -£.iL)- 

2 b 2 b 2 M 2 2 


~ tan $ = 0 =e 

b 


When $ is specified, this equation is solved so that it satis- 
fies e < 10 " 4 . The corresponding Mach number component Mjj can 

then be solved with equation (36). 
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•Hie velocity vector the medium range in the aerodynamic 
coordinate system is given as follows in detail according to 
equation (45): 


Ji 

|U« 


i [cos a cos y (AUBI)tcos a sin y ( g C M 4 e ) cos ip) 

+ sin a { g (M 4 e) sin ip)] 

♦j [cos y (g(M; e ) cos <p)-sin y (A1 + B1 )] 

♦k [sin a cos y (AW Bl)*sin a sin y(g( M 4e ) cos <p)-cos a cos <pj 


A 1 and B 1 are given by 


A1 = cos a cos y (a -b tan y ) 
B1=sin a cos v (c -b ■— ”) 


(45a) 


APPENDIX V 

In the case where the internal shock is to produce an addi- 
tional component in the direction of the leading edge, then 
U 3 ojit 0=»K. • Then the velocity in the medium range is 

given by equation (45a), and in addition to A 1 and B 1 we have 

Cl s (b-b*K) cos v 

Therefore, the sum (A 1 + B l) in equation (45a) is the expanded 
into (A 1 + B 1 -C 1). 
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TABLE la: THEORETICAL HANDLING OF LEE SIDE FLOW. 
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TABLE lb: THEORETICAL TREATMENT OP LEE SIDE PLOW 
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Parabolic lino Expansion Reflected expansion waves Mach cone 




















TABLE 2: SUMMARY OF MOST IMPORTANT EXPERIMENTS ON LEE SIDE FLOW 
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Fig. 


Fig. 


Subsonic loading edge 



Primary vortex 

L_ 1 .fy 

I Secondary vortex 

i I 


Supersonic leading edge 


l u - iu- 

V * 



a) attached b) Leading c) attached d) shock-induced 

flow edge sepa- flow separation 

ration 


1. Plow models on the lee side of slender delta wings. 



Mach number perpendicular to 
leading edge 

2. Flow regions over the thin flat delta wing. 
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SA St, 



Section perpendicular to leading edge (VK 

A - attached flow 

SL- separation line 

SA- secondary separation 


Vortex formation and 


stream line course 



Pressure distribution 

Pressure part of secondary vortex 


Fig. 3. Lee side flow over slender delta wing with subsonic 
leading edge. 






Pressure distribution 



Pig. 4. Lee side flow over slender delta wing with supersonic 
leading edge. 
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Attached boundary layer' ^ 

4 with flow perpendicular^ wonlcai flow 
to central plane ' 
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Attached leading edge shock") 

Internal shock ( 

Prandtl-Meyer expansion rNon-conlcal flow 


Vortex 


^ Attached flow \ 

a / rr — Internal shock . J* 

A /, « Prandtl-Meyer expansion 
Vortex -J 


Non-conlcal flow 



Fig. 5. The hypersonic flow inode! of Whltehead/Heffner/Rao [6]. 




Attached flow ) Conical flow 

.interference of leading edge 
shock and separation bubble leads 
■ to leading edge separation * 

'Separated leading edge shock 


Non-coni cal flow 


Attached shock 
Free vortices 
Internal Shock 
/ Separation bubble with 
constant width J 

Separated leading edge shock 
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Fig. 6. The supersonic flow model of Rein [46] for separated 
leading edge shock. 
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Mach cone 


.Internal shocks 


Leading edge shock 


Boundary layer 

'Separation 


Mach cone 
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Leading edge shock, 
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Fig. 7. Flow models for delta wings in the hypersonic regions 
for medium angles of attack. 
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Mach number perpendicular to leading edge Smalls 


Secondary vortex 

Underside flow with 
1 stagnation line 



Leading edge vortex 
x Separated leading edge shock 


Fig, 8. The transition from shock-induced separation to leading 
edge vortex formation. 
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Slenderness s/1 * 0.31 

1 Convex delta model with 
leading edge angle tp ■ 55° 

2 Flat delta model with 
leading edge angle ^ “ 40° 

3 Flat cone model with 
leading edge angle tp ** 75 



Fig. 10. Wind tunnel Model I 
(CUED, RAE) 


Slenderness s/1 * 0.30 

4 Flat delta model with 
leading edge angle tp * 38.9 

5 Flat cone model with 
leading edge angle ip 85 75 

6 Nonweiler Wave rider with 
leading edge angle ip * 39.5° 
(initial parameter M 0 -« 


Fig. 11. Wind tunnel Model II 
(ILR) 
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Slenderness s/1 =0.3 



Fig. 12. Modified cone model (ILR) 


Gamb ridge wind 




ILR wind 
tunnel 



Fig. 13. 


Model suspension in 
various wind tunnels 


^ ~ £ 


• M.f 1.01 



Fig. 14. Pitot probes, cali- 
bration as a function 
of Incident flaw angle. 


Angle of attack 0 



e, = position of attachd flow 

0p: = primary vortex center 

°« = secondary vortex center Fig . 15< De f in ition of oil film 

= primary vortex intensity i TOage evaluations 

♦ s « secondary vortex intensity 
= primary vortex parameter 
= secondary vortex parameter 
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*: Separation with leading edge vortex “ 

» O 



Fig. 


16 .Angle of attack and Mach number 
perDendicular to leading edge 
for delta wings with different 
cross-section shapes. 


nic 




Fig. 18. Isobars on deltas 
shaped lee side of 
the Model I, 


M - 3.5 
oo 





Fig. 19. Flow field analysis over the flat lee side of the delta 
wing at a = 10° f M eo = 2.5 
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Fig. 22. Position of primary CPZ) 
and secondairy (SZ) vor- 
tex overvthe delta wing. 
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Fig. 23. Pressure coefficients 
over span for differ- 
ent Mach numbers, 

a- 10° 

( a N =20°) 




i 





Fig.. 24. Changes of minimum and maximum pressure on lee side of delta wings 

as a function of Mach number, , A o 

/ <X= IQ 



Fig. 25. Oil film and Schlieren photographs with the flat lee side for 
different Mach numbers, ^ _ jqO 


■ * ... 




Fig. 27 . Vortex positions when crossing the Stanhrook-Squire region 



Oil film 


Fig. 28, Flow field analy- 
sis over flat lee 
side of the delta 
wing at 


«= 10° , M*. =* 3.5 
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Fig. 29, Experimentally determined flow model, M 
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Fig. 3Q. Oil film photographs 
on the flat lee side 
of the delta wing 
for different angles 
of attack, M “3.5 
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Fig. 31. 
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Pressure coefficients over 

span for different angles 

of attack, M * 3.5. 
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Fig. 32. Vortex position for different 
angles of attack 


Angle of attack 13 



Fig. 


33. Vortex intensity repre- 
sented by delta flow 
angle in the oil film 
photograph for different 
Mach numbers 
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Fig. 35. Collection of experi 
ments on the delta- 
shaped lee side 

M oo = 2 * 5, a= J0° 
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Pitot Isobars on the top side of a delta wing for different angles 

of attack. Measurement plane */l =0.8, M =2.6 (M. = 1.39) 
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Relative half span y/s 


Fig. 37. Pressure coefficients over span for different Mach numbers 

a =10° ( a N = 17°) 


Conical flow 


Transition region 
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Fig. 38. Model of the flow over the delta— shaped lee side for the example 

of M = 3.5 and a = 10° 
oo 
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Fig. 39. Pressure coefficients- over span for different angles of attack 

M »3.5 
00 



Position of vortex center 
as a function of angle of 
attack. 


138 



a.. K>* 


II* 


II* II* 



Relative half span y/s 

Tig. 41. Oil film photographs and pressure distributions with different 
top side shapes and fixed incident flow conditions 

=3.5 , a~10° 
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Fig. 42. Position of vortex center and attachment line as a function of 
shape parameter 0 ' 
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Fig. 43: Influence of 

extreme leading 
edge angles for 
the example of 
the cone model, 


CL= 12° 


'll 

Shock, (separation) 
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Or' G/ 4/4, 
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Fig. 44. Flow around the modified cone model shown with pressure distributions 
and flow models-, M = 3.06. 

7 m 
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Fig. 45. 


Comparison of flow 
fields of two delta 
wings with Identical 
lee side but differ- 
ent bottom side 
shape * 2.6. 



Fig. 46. 



Position of stagna- 
tion lines on the 
flat pressure side 
of the delta wing 
in the region of 
separation with in- 
ternal shock. 
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Fig. 47. Vortex inten- 
sity as a func- 
tion of angle 
of attack, 

M * 3.5. 
oo 




Pig. 48. Comparison of pressure distributions on the lee side with different 
cross-section shapes, M =3.5. 
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Fig. 50 
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delta wings, *1 = 2,6 
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Fig. 51. Restrictions for the design of the top side of a delta wing for 
cruise flight. 
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Fig. 53. Geometry of the simplified lee 
side flow model. 
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Fig. 54. Flow Model 


Base point of internal shock 
Leading edge shock 


Section A-A 


145 



Angle of attack and model values 
i M..0, I 


Leading edge 
type 


Supersonic L.E. 


Subsonic - L.E. 




Attached shock 


Leading edge' 
shock 


separated shock 


Flow perpendicular 
to leading edge 


Approx, solution for 
— boundary ‘layer L — 

I GS.f(M.,aJ 


Approx, solution for 
shock position and shape 

Oblique shock , 

relationships „ 


_Flow perpendicular 
to leading edge 


j frandlUM*y«r j 

expansion over B.L. 


First characteristic, 
which has to intersect 
the Internal shock 


Output j 


Initial point of 
L internal shock 

f— 

Iteration of 
internal shock J 


-Influencing of shock 
I* - Intensity. 1 


Supersonic 


j Frond [ 

expansion over B.E. 


Boundary layer thickening 
behind Internal shock 


Mach number component 
''behind internal s"hock. 



JL 

Flow components 
'perpendicular to* 
symmetry line"^ 


Subsonic 


Panel method 


V; * 0 


Fig. 55. Macro-flow diagram for calculating lee side flow. 
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Fig. 56. Coordinates and unit vector system. 



Thickness parameter C 


Fig. 57. Flow region represented in the C= f (SI) diagram for 
the "flat delta models." 
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Fig. 58. Leading edge shock 

profiles for different 
angles of attack, 


? o.. 


"Ov , Mach cone 



^Estimated variation of 
\ leading edge shock 


Fig. 59. 


Profile of leading edge 
shock compared with 
different theories, 

M «= 3.5 
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Moeckel [83] 
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Fig. 60. Separation of the separated 
leading edge shock as a 
function of Mach number. 
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Fig. 61. Angle of leading edge 
shock as a function of 
Mach number. 
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Fig. 62. Geometry for the separated 
leading edge shock. 
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Pig. 64. Solution of the displacement problem using the Panel method. 
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Fig. 65. Influence of boundary layer displacement thickness on pressure 
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Fig. 66: Influence of inclination angle v s of leading edge shock on 

pressure coefficients and position of internal shock, “ 3.55, 

a = 14.5° 
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Fig. 67: Influence of trace angle S of the leading edge shock on pressure 

coefficients and position of internal shock, = 3.55, a = ]4.5° 
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Fig. 72: Comparison theory - 

experiment: Pitot Iso' 
bars measurement of 
Monnerie and Werla, 
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Fig. 73. Influence of Mach number and angle of attack variation on the 
flow field of the lee side of delta wings with A. = 72.5° 
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